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FORTRAN PROGRAM FOR CALCULATING VELOCITIES AND STREAMLINES 
ON THE HUB-SHROUD MID-CHANNEL FLOW SURFACE OF AN 
AXIAL- OR MIXED-FLOW TURBOMACHINE 
I - USER'S MANUAL 

by Theodore Katsanis and William D. McNally 
Lewis Research Center 

SUMMARY 

A FORTRAN-IV computer program has been developed which obtains a subsonic or 
transonic, nonviscous flow solution on the hub-shroud mid-channel flow surface of a 
turbomachine. The flow must be essentially subsonic, but there may be locally super- 
sonic flow. The solution is for two-dimensional, adiabatic shock-free flow. The blade 
row may be fixed or rotating, and may be twisted and leaned. The flow may be axial or 
mixed, up to approximately 45° from axial. Upstream and downstream flow conditions 
can vary from hub to shroud, and provision is made for an approximate correction for 
loss of stagnation pressure. 

The basic analysis is based on the stream function and consists of the solution of 
the simultaneous, nonlinear, finite-difference equations of the stream function. This 
basic solution, however, is limited to strictly subsonic flow. When there is locally 
supersonic flow, a transonic solution must be obtained. The transonic solution is ob- 
tained by a combination of a finite-difference stream -function solution and a velocity- 
gradient solution. The finite -difference solution at a reduced mass flow provides infor- 
mation which is used to obtain a velocity-gradient solution at the full mass flow. 

The program input consists of blade and flow -channel geometry, upstream and 
downstream flow conditions from hub to shroud, and mass flow. The output includes 
streamline coordinates, flow angles, and velocities on the mid-channel flow surface; in- 
cidence and deviation angles at the blade leading and trailing edges; and approximations 
to the blade surface velocities. The output may also include input information for a 
blade -to -blade flow analysis program. 

The program is reported in two volumes, with part I as the user's manual and 
part II as the programmer's manual. This report, part I, contains all the information 
necessary to use the program as is. It explains the equations involved and the method of 
solution and gives a numerical example to illustrate the use of the program. Part II in- 
cludes the complete program listing and a detailed program procedure. 


INTRODUCTION 


The design of blades for compressors and turbines ideally requires analysis methods 
for unsteady, rotational, three-dimensional, viscous flow through a turbomachine. 
Clearly, such solutions are impossible at the present time, even on the largest and fast- 
est computers. The usual approach at present is to analyze only steady flows and to 
separate inviscid solutions from viscous solutions. Three-dimensional inviscid solutions 
are just beginning to be contemplated for coming generations of computers. So at pres- 
ent, inviscid analyses usually involve a combination of several two-dimensional solutions 
on intersecting families of stream surfaces to obtain what is called a quasi-three- 
dimensional solution. 

Since there are several choices of two-dimensional surfaces to analyze, and many 
ways of combining them, there are many approaches to obtaining a quasi-three- 
dimensional solution. Most two-dimensional solutions are either on a blade-to-blade 
surface of revolution (Wu's Sj surface, ref. 1) or on the meridional or mid-channel 
stream surface between two blades (Wu's S£ surface). However, when three- 
dimensional effects are most important, significant information can often be obtained 
from a solution on a passage cross-sectional surface (normal to the flow). This is 
called a channel solution (see fig. 1). 



Figure 1. - Two-dimensional analysis surfaces in a turbomachine. 
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In this report a solution to the equations of flow on the meridional S2 surface is 
carried out. This solution surface is chosen when the turbomachine under consideration 
has significant variation in flow properties in the hub-shroud direction. A solution on 
the meridional surface will show this variation. The solution can be obtained either by 
the quasi-orthogonal method, which solves the velocity-gradient equation from hub to 
shroud on the meridional flow plane (ref. 2), or by a finite-difference method, which 
solves a finite-difference equation for stream function on the same flow plane. The 
quasi-orthogonal method is efficient in many cases and can obtain solutions into the tran- 
sonic regime. However, there is difficulty in obtaining a solution when aspect ratios 
are above 1. Difficulties are also encountered with curved passages and Iow-hub-tip- 
ratio blades. For such cases, the most promising method is the finite-difference solu- 
tion, but this solution is limited to completely subsonic flows. 

Two finite-difference programs for flow on the mid-channel surface of a turbo- 
machine have been reported in the literature (refs. 3 and 4). Since both are finite- 
difference methods, they are necessarily limited to subsonic flow cases. Marsh's 
method (ref. 3), termed the matrix throughflow method, closely follows the development 
given by Wu in reference 1. However, the computer program was not included in refer- 
ence 3, nor is it available to the general user. Davis' program is provided in refer- 
ence 4 but is limited to certain families of compressor blades and flow surfaces. 

The method described in this report uses both the finite-difference and the quasi- 
orthogonal (velocity gradient) methods, combined in a way which takes maximum advan- 
tage of both. The finite-difference method is used to obtain a subsonic-flow solution; 

The velocity-gradient method is then used if necessary to extend the range of solutions 
into the transonic regime. 

A computer program, MERIDL, has been written to perform these calculations. 

This program is written for axial- or mixed-flow turbomachines, both compressors and 
turbines, to approximately 45° from axial. Upstream and downstream flow conditions 
can vary from hub to shround. The solution is for compressible, shock-free flow, or 
incompressible flow. Provision is made for an approximate correction for loss of stag- 
nation pressure through the blade row. The blade row may be either fixed or rotating 
and may be twisted and leaned. The blades can have high aspect ratio and arbitrary 
thickness distribution. 

The solution obtained by this program also provides the information necessary for • 
a more detailed blade shape analysis on blade-to-blade surfaces (fig. 1). A useful pro- 
gram for this purpose is TSONIC (ref. 5). Information needed to prepare all the input 
for TSONIC is calculated and printed by MERIDL. 

The MERIDL program has been implemented on the NASA Lewis time- sharing 
IBM-TSS/360-67 computer. For the numerical example of this report, storage of vari- 
ables required 60 000 words for a 21 x 41 grid of 861 points. Variable storage could be 
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easily reduced by equivalencing of variables or by using a coarser mesh. Storage for 
the program code is 18 000 words. This storage could be reduced by overlay of code. 
Run times for the program range from 3 to 15 minutes on IBM 360-67 equipment, de- 
pending upon the mesh size used and the compressibility of the flow. 

The MERIDL program is reported in two volumes, with part I as the user's manual 
and part II (ref. 6) as the programmer's manual. This report, part I, contains all the 
information necessary to use the program as is. It explains the method of solution, de- 
scribes the input and output, gives a numerical example to illustrate the use of the pro- 
gram, and derives the equations used (in the appendixes). Part II includes a complete 
program listing, detailed program procedure, and appendixes which derive special nu- 
merical techniques used. 


METHOD OF ANALYSIS 
BASIC ASSUMPTIONS 

It is desired to determine the flow distribution through a stationary or rotating cas- 
cade of blades on a mid-channel hub-shroud stream surface. The following simplifying 
assumptions are used in deriving the equations and in obtaining a solution: 

(1) The flow relative to the blade is steady. 

(2) The fluid is a perfect gas with constant specific heat C p . 

(3) The fluid is a nonviscous gas. 

(4) There is no heat transfer. 

(5) The mid-channel surface is a stream surface which has the same shape as the 
blade mean camber surface, except near the leading and trailing edges, where an arbi- 
trary correction is made to match the free-stream flow. 

(6) The only forces are those due to momentum and pressure gradient. 

(7) The velocity varies linearly between blade surfaces. 

(8) The relative stagnation pressure loss is known through the blade row. 

The flow may be axial or mixed, to about 45° from axial. There may be a variation 
of whirl, stagnation pressure, and stagnation temperature from hub to shroud, both up- 
stream and downstream of the blade row. The blade row may be either fixed or rotating, 
with leaned and twisted blades. Within the given assumptions, no terms are omitted 
from the equations. 
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SOLUTION BY COMBINATION OF METHODS 

A flow analysis on the meridional flow surface can be obtained either by the velocity- 
gradient method or by the finite-difference method. The finite-difference method is 
limited to subsonic flow, whereas the velocity-gradient method is limited to relatively 
low-aspect-ratio blades. The most accurate solution is obtained by the finite-difference 
technique, so that this method is used where possible (i. e. , for subsonic flow). With 
locally supersonic flow, the finite-difference solution is first obtained at a reduced mass 
flow for which the flow field is completely subsonic. The streamline curvatures and flow 
angles throughout the passage which are obtained from this solution provide the informa- 
tion necessary to obtain an approximate velocity-gradient solution at full mass flow, re- 
gardless of aspect ratio. 


SUBSONIC STREAM-FUNCTION SOLUTION 


The stream -function equation is a partial differential equation on a mid-channel hub- 
shroud stream surface (see assumption 5). This equation is in one unknown (the stream 
function) as a function of two variables, r and z (see fig. 2). This is equation (Al) and 



Figure 2. - Cylindrical coordinate system and velocity components. 
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is derived from Wu's equation (eq. (107a), ref. 1) for the stream function on what he 
calls an S 2 surface. Equation (Al) is nonlinear but can be solved iteratively by the 
finite-difference method when the flow is completely subsonic. 

A finite region (as indicated in fig. 3) is considered for the solution of equation (Ai). 
It is assumed that the upstream and downstream boundaries are sufficiently far from 
the blade so as to have a negligible effect on the solution. Equation (Al) is elliptic for 
subsonic flow. Therefore, when the flow is entirely subsonic, equation (Al) can be 
solved when proper boundary conditions are specified on the entire boundary of the re- 
gion. These conditions are the values of the stream function on all four boundaries. 

The stream function has the value 0 on the hub and 1 at the shroud. The value of the 
stream function on the upstream and downstream boundaries can be calculated if the 
stagnation pressure, stagnation temperature, and whirl distribution from hub to shroud 
are specified upstream and downstream of the blade. 

The numerical solution of equation (Al) is obtained by the finite-difference method. 

A grid must be used for the finite-difference equations. The type of grid that was 
chosen is an orthogonal mesh which is generated by the program. The method used to 
generate the grid follows that reported in reference 7. The space between the hub and 
the shroud is divided into equal increments along several radial lines. Spline curves 
are then fit through the resulting points to obtain the streamwise orthogonals (see fig. 4). 



The normal orthogonals are obtained by a predictor-corrector technique. This tech- 
nique is analogous to the second-order Runge-Kutta method for solving ordinary differ- 
ential equations, known also as the improved Euler method or Heun’s method (ref. 8). 
The orthogonal mesh coordinates are s and t. The s-coordinate is in the streamwise 
direction and the t-coordinate is normal to this, as indicated in figure 4. With the 
mesh determined, the finite-difference equations can be derived. The finite-difference 
equations on the orthogonal mesh are given in part II of this report. 

The finite-difference equations are nonlinear since the original equation (Al) is 
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nonlinear. These equations can be solved iteratively. On the first iteration an initial 
density is assumed; this linearizes some of the terms. The remaining nonlinear terms 
are omitted for the first iteration so that the finite-difference equations are entirely 
linearized. These linearized equations are then solved to obtain the first approximate 
solution for stream function. This solution provides information used to obtain a better 
estimate of the density and an estimate of the other nonlinear terms. The equations 
are then solved again to obtain an improved solution. This process is repeated, and 
by iteration a final converged solution can be obtained if the flow is subsonic. 

For each step of this iteration the linearized finite- difference equations must be 
solved. The method used to solve the equations is successive overrelaxation (ref. 9) 
with an optimum overrelaxation factor. Since this is also an iterative method, we have 
two levels of iteration. The overrelaxation is performed in the "inner iteration, " and 
the corrections to the nonlinear terms are made in the "outer iteration. " 

After the stream function is obtained, the velocity distribution is obtained by nu- 
merical partial differentiation of the stream function and by using equations (A5) 
and (A6). The details of the numerical procedure and programming technique are de- 
scribed in part II. 


TRANSONIC VELOCITY-GRADIENT APPROXIMATE SOLUTION 

For the case where there is locally supersonic flow, equation (Al) is no longer el- 
liptic in the entire region but is hyperbolic in the region of supersonic flow (ref. 10). 
This changes the boundary conditions and means that there will probably be shock losses 
in going from supersonic to subsonic flow. The finite-difference method cannot be used 
with locally supersonic flow. However, an approximate solution can be obtained by 
getting a reduced-flow solution with the finite-difference method and extending this to the 
full flow by using the velocity-gradient method. This technique is described in refer- 
ence 5. 

The velocity-gradient equations are equations (A7) to (All). Equation (A7) is solved 
as an initial-value problem, where the velocity W is specified at the hub for any given 
vertical mesh line running from hub to tip. By finding several solutions for varying 
values of W at the hub, a solution satisfying the specified mass flow (eq. (A12)) will be 
found. When equation (Al) has been solved, subject to giving the correct mass flow, for 
every hub- shroud mesh line in the region, the entire velocity distribution at full mass 
flow has been obtained. 
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BLADE SURFACE VELOCITIES 


The solution which is obtained by either the finite-difference or velocity -gradient 
method is for the mid-channel surface between the blades. Of greater interest are the 
blade surface velocities. These can be estimated since the blade loading is dependent 
on the rate of change of whirl. By assuming a linear variation of velocity between blade 
surfaces, equation (A13) can be derived for calculating the blade surface velocities. 


APPLICATION OF PROGRAM 


The program can be used both for analysis and as a design tool. When used for de- 
sign, other programs should be used with this program. For axial compressors, refer- 
ence 11 describes a program that will give blade me an- camber -line coordinates and 
thicknesses for an axial compressor blade. This blade design can be checked by using 
the MERIDL program to analyze the flow distribution in detail. Usually, changes must 
be made to the blade design to achieve a desirable flow distribution. These changes may 
involve more than just the blade shape; for example, hub and shroud profile, inlet and 
outlet whirl distribution, and loss distribution may have to be changed. Of course, the 
accuracy of the MERIDL solution depends on the accuracy of the boundary conditions 
used. 

When a reasonable flow pattern is achieved by the MERIDL flow analysis on the mid- 
channel flow plane, more detailed blade surface velocities can be obtained on blade-to- 
blade planes by flow analyses on various blade elements from hub to shroud. A useful 
program for this purpose is TSONIC (ref. 5). Most of the information required to com- 
pute input for TSONIC is calculated and printed directly by MERIDL. The section 
Printed Output explains how to obtain and use this output from MERIDL. Further 
changes in blade shape or whirl distribution may be considered at this time. Refer- 
ence 12 (Ch. VII) gives information on incidence and deviation for good design and for 
off-design conditions. 

For cases when the flow is well guided in the channel but has large variations, both 
blade to blade and hub to shroud, the CHANEL program (ref. 13) is useful. The 
CHANEL program obtains a solution on a channel cross-section surface. The CHANEL 
•program is particularly useful for calculating choking mass flow through a blade row. 


DESCRIPTION OF INPUT AND OUTPUT 


The principal block of input required by the program is a geometrical description of 
the blade row to be analyzed, given in the form of blade sections from hub to shroud. 
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Each blade section is described by a set of z-, r-, 0-, and normal-thickness coordi- 
nates on any general, smooth surface of revolution. There does not have to be any 
geometrical relation between analogous points on adjacent blade sections. Other inputs 
include upstream and downstream flow conditions, a geometric description of hub and 
shroud, appropriate gas constants, operating conditions such as mass flow and rotational 
speed, and a description of the finite -difference solution mesh. Requests are also given 
for various blocks of output data. 

Output is given at any or all of three principal geometric locations: (1) on all mesh 
points of the orthogonal solution mesh, (2) along user-designated streamlines through 
the blade row, or (3) along straight lines from hub to shroud (station lines) located either 
inside or outside the blades or at the blade leading or trailing edges. In each of these 
locations, output consists mainly of z- and r-coordinates of orthogonal mesh lines or 
streamlines, stream function, relative velocities and velocity components, flow angles, 
and streamline curvature. Along streamlines and station lines which lie within the 
blades, estimates are also given of the blade surface velocities. Other desired output 
may be easily obtained by user modification of the program. 


INPUT 

Figure 5 shows the input variables as they are punched on the data cards. The first 
input data card is for a title, which serves for problem identification. Any information 
may be put in the 80 columns of this card. 

All the numbers on the three input cards beginning with MBI, LSFR, and IMESH are 
integers (no decimal point) in a five-column field (see fig. 5). These must all be right- 
adjusted. The input variables on all other data cards are real numbers (punch decimal 
point) in 10-column fields. 

Figure 5 indicates that several options exist for the statement of upstream and 
downstream flow conditions. First, the user can specify either whirl or absolute tangen- 
tial velocity (LAMIN or VTHIN, and LAMOUT or VTHOUT). Whirl must be given as a 
function of stream function (SFIN and SFOUT). Tangential velocity is usually given as a 
function of radius (RADIN and RADOUT). And finally, on the downstream boundary, 
either absolute total pressure or absolute total pressure loss may be given (PROP or 
LOSOUT). 

Input variables are both geometric and nongeometric. The geometric input variables 
are shown in figures 6 to 10. Further information concerning the input variables is 
given in the section Special Instructions for Preparing Input, pages 19 to 24. 
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Figure 5. - Input form. 
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ZTOUT 


r 



Figure 7. - Input variables - 


□ Hub or shroud spline point 
0 Hub point where orthogonal 
mesh begins, ends, or 
changes spacing 



■ hub, shroud, and blade sections. 
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0 Hub point where orthogonal mesh begins, 
ends, or changes spacing 



ZOMBO 


ZOMOUT 




e 



Figure 10. - Input variables - blade section. RBL must be given for each ZBL.THBL 
location (see fig. 6). 


Input Dictionary 


The input variables are described in terms of a consistent set. of metric units: new- 
tons, kilograms, meters, joules, kelvins, and seconds. The program, however, will 
run with input in any consistent set of units. 

The input variables are the following: 


GAM 

AR 

MSFL 

OMEGA 

RED FAC 


VELTOL 


Specific -heat ratio, y 

Gas constant, J/(kg)(K) 

Total mass flow through entire circumferential annulus of machine, kg/sec 

Rotational speed, w, rad/sec. Note that u> is negative if rotation is in op- 
posite direction of that shown in fig. 2. 

Factor by which mass flow (MSFL) must be reduced in order to assure sub- 
sonic flow throughout flow passage. REDFAC may be left blank, in which 
case a value of 1.0 will be used. See section (k), p. 24. 

Convergence tolerance on maximum velocity change in each outer iteration, 
over all mesh points, for reduced mass flow. VELTOL may be left blank, 
and a value of 0. 01 will be used by the program. Whatever value is given 
is multiplied by the minimum of FNEW and DNEW before it is printed and 
subsequently used by the program. A value of 0. 001 for VELTOL is a 
tight tolerance, 0. 01 is a medium tolerance, and 0. 1 is a loose tolerance. 
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FNEW 

DNEW 

MBI 

MBO 

MM 

MHT 

NBL 

NHUB 

NTIP 

NIN 

NOUT 

NBLPL 

NPPP 


Damping factor on calculation of . from outer iteration to outer iteration. 

A value of 0. 5 is suggested for FNEW. FNEW may be left blank, in which 
case the program will use a value of 0. 5. See section (b), p. 20. 

Damping factor on calculation of 3(rVg)/0r within blade row from iteration 
to iteration. DNEW may be left blank, and the program will use a value of 
0. 5. DNEW does not have to be equal to FNEW. See section (b), p. 20. 

Number of vertical mesh lines from left boundary of orthogonal mesh 
(ZOMIN) to point of first mesh size change (ZOMBI). See fig. 8 and sec- 
tion (d), p. 21. 

Total number of vertical mesh lines from left boundary of orthogonal mesh 
(ZOMIN) to point of second mesh size change (ZOMBO). See fig. 8 and 
section (d), p. 21. 

Total number of vertical mesh lines from left to right boundaries of orthog- 
onal mesh (ZOMIN to ZOMOUT), maximum of 100. See fig. 8 and sec- 
tion (d), p. 21. 

Total number of horizontal mesh spaces from hub to shroud of orthogonal 
mesh, maximum of 100. See fig. 8 and section (d), p. 21. 

Number of blades in total circumference of blade row. 

Number of spline points given in ZHUB and RHUB arrays, maximum of 50. 

See fig. 6 and section- (c), p. 20. 

Number of spline points given in ZTIP and RTIP arrays, maximum of 50. 

See fig. 6 and section (c), p. 20. 

Number of data points given in upstream arrays of flow properties (SFIN, 
RADIN, TIP, PRIP, LAMIN, VTHIN), maximum of 50. See fig. 7 and 
section (e), p. 21. 

Number of data points given in downstream arrays of flow properties (SFOUT, 
RADOUT, PROP, LOSOUT, LAMOUT, VTHOUT), maximum of 50. See 
fig. 7 and section (e), p. 21. 

Number of blade planes or blade sections on which data (ZBL, RBL, THBL, 
TNBL) are given to describe mean flow surface and blade thickness, max- 
imum of 50. See fig. 6 and section (f), p. 22. 

Number of data points per blade section or blade plane in ZBL, RBL, THBL, 
and TNBL arrays, maximum of 50. See fig. 10 and section (f), p. 22. 
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NOSTAT 


NSL 


LSFR 

LTPL 

LAMVT 

ZOMIN 

ZOMBI 

ZOMBO 

ZOMOUT 

ZHUB 

RHUB 

ZTIP 

RTIP 


Number of hub-shroud stations (located by coordinates in ZHST and ZTST) 
at which output is desired, maximum of 50. See fig. 9 and section (h), 
p. 23. NOSTAT may be left blank, in which case no input cards should be 
included for ZHST and ZTST arrays. 

Number of streamlines from hub to shroud (designated by values in FLFR) 
at which output is desired, maximum of 50. See fig. 9. NSL may be left 
blank, in which case no cards should be included for FLFR array. If NSL 
is left blank, the program will set it equal to 11 and print requested stream- 
line output on 11 streamlines which vary by 10 percent of total flow (i. e. , 

0, 10, 20, . . . , 100 percent). 

Integer (0 or 1) indicating whether upstream and downstream flow conditions 
are given as a function of stream function (0) or radius (1). If 
LAMVT = 0, LSFR must equal 0. 

Integer (0 or 1) indicating whether downstream total pressure (0) or frac- 
tional loss of stagnation pressure (1) is given as input. 

Integer (0 or 1) indicating whether upstream and downstream whirl (0) or 
tangential velocity (1) is given as input. 

z -Coordinate of intersection of left boundary of orthogonal mesh with hub 
profile, m. See figs. 6 and 8 and section (d), p. 21. 

z-Coordinate of intersection of vertical mesh line with hub profile where 
first change in mesh spacing occurs (MBI), m. See figs. 6 and 8 and 
section (d), p. 21. 

z-Coordinate of intersection of vertical mesh line with hub profile where 
second change in mesh spacing occurs (MBO), m. See figs. 6 and 8 and 
section (d), p. 21. 

z-Coordinate of intersection of right boundary of orthogonal mesh (MM) with 
hub profile, m. See figs. 6 and 8 and section (d), p. 21. 

Array of z-coordinates of input points defining hub or bottom boundary of 
flow channel, m. See fig. 6 and section (c), p. 20. 

Array of r-coordinates of input points defining hub or bottom boundary of 
flow channel, m. See fig. 6 and section (c), p. 20. 

Array of z-coordinates of input points defining shroud or top boundary of 
flow channel, m. See fig. 6 and section (c), p. 20. 

Array of r-coordinates of input points defining shroud or top boundary of 
flow channel, m. See fig. 6 and section (c), p. 20. 
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ZfflN 

ZTIN 

SFIN 

RADIN 

TIP 

PRIP 

LAMIN 

VTHIN 

ZHOUT 

ZTOUT 

SFOUT 

RADOUT 

PROP 


z -Coordinate of intersection with hub profile of line on which upstream flow 
conditions are given, m. See fig. 7 and section (e), p. 21. 

z -Coordinate of intersection with shroud profile of line on which upstream 
flow conditions are given, m. See fig. 7 and section (e), p. 21. 

Array of values of stream function for input points from hub to shroud along 
line on which upstream flow conditions are given. See fig. 7 and sec- 
tion (e), p. 21. 

Array of r-coordinates of input points along line from hub to shroud on which 
upstream flow conditions are given, m. See fig. 7 and section (e), p. 21. 

Array of absolute total temperatures TJ at input points along line from hub 
to shroud on which upstream flow conditions are given, K. See fig. 7 and 
section (e), p. 21. 

Array of absolute total pressures p! at input points along line from hub 

1 2 

to shroud on which upstream flow conditions are given, N/m . See fig. 7 
and section (e), p. 21. 

Array of values of absolute whirl frVgj at input points along line from hub 

i 2 

to shroud on which upstream flow conditions are given, m /sec. See 
fig. 7 and section (e), p. 21. 

Array of values of absolute tangential velocity at input points along 

line from hub to shroud on which upstream flow conditions are given, 
m/sec. See fig. 7 and section (e), p. 21. 

z-Coordinate of intersection with hub profile of line on which downstream 
flow conditions are given, m. See fig. '7 and section (e), p. 21. 

z-Coordinate of intersection with shroud profile of line on which downstream 
flow conditions are given, m. See fig. 7 and section (e), p. 21. 

Array of values of stream function for input points from hub to shroud along 
line on which downstream flow conditions are given. See fig. 7 and sec- 
tion (e), p. 21. 

Array of r-coordinates of input points along line from hub to shroud on 
which downstream flow conditions are given, m. See fig. 7 and sec- 
tion (e), p. 21. 

Array of absolute total pressures p' at input points along line from hub to 

o 2 

shroud on which downstream flow conditions are given, N/m . See fig. 7 
and section (e), p. 21. 
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LOSOUT 

LAMOUT 

VTHOUT 

ZBL 

RBL 

THBL 

TNBL 


ZHST 

ZTST 


Array of fraction of absolute total pressure loss (p^ - p’ 0 )/p’ 0 at input 

points along line from hub to shroud on which downstream flow conditions 
are given. See fig. 7 and section (e), p. 21. 

Array of values of absolute whirl (rVgJ at input points along line from hub 

o 2 

to shroud on which downstream flow conditions are given, m /sec. See 
fig. 7 and section (e), p. 21. 

Array of values of absolute tangential velocity (VgJ at input points along 

' o 

line from hub to shroud on which downstream flow conditions are given, 
m/sec. See fig. 7 and section (e), p. 21. 

Two-dimensional array of z-coordinates of points describing mean blade sur- 
face, m. See fig. 10 and section (f), p. 22. This surface is described by 
a series (from 2 to 50) of blade sections from hub to shroud. The hub sec- 
tion is given first, followed by successive sections up to the shroud. 

Two-dimensional array of r-coordinates, corresponding to ZBL, of points 
describing mean blade surface, m. See fig. 6 and section (f), p. 22. 

Two-dimensional array of 0-coordinates, corresponding to ZBL, of points 
describing mean blade surface, rad. See fig. 10 and section (f), p. 22. 
Theta is positive in direction of positive rotation (see fig. 2). The 
origin of 0-coordinates can be anywhere around the circumference. 

Two-dimensional array of blade normal thicknesses, corresponding to ZBL, 
RBL coordinates, m. See fig. 10 and section (f), p. 22. TNBL is 
thickness on a surface of revolution and is normal to blade mean camber 
line. It is not the thickness at a constant value of z-coordinate from suc- 
tion surface to pressure surface (see fig. 10). Thickness from either a 
conical blade section or a cylindrical blade section can be used for TNBL; 
it makes little difference. This input is used only to calculate local blade 
blockage and usually has a minor effect on computed results. 

Array of z-coordinates of intersections of hub-shroud output station lines 
with hub profile, m. See fig. 9 and section (h), p. 23. No input cards 
should be given for ZHST and ZTST if user does not wish output on hub- 
shroud station lines (i. e. , NOSTAT = 0). 

Array of z-coordinates of intersections of hub-shroud output station lines 
with shroud profile, m. See fig. 9 and section (h), p. 23. 
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FLFR 


Array of values of stream function designating stream lines along which out- 
put is to be printed. See fig. 9 and section (h), p. 23. If no cards are 
given for FLFR array (i.e. , NSL = 0), the program will automatically 
assign the following 11 values to FLFR: 0.0, 0.1, 0.2, ..., 0.9, 1.0. 

The remaining seven integer variables, beginning with IMESH, are used to indicate 
what output is desired. The program reduces the input mass flow by REDFAC and 
solves the resulting problem iteratively with a stream -function analysis. The results of 
this analysis are then used with full mass flow to obtain a transonic analysis by velocity- 
gradient methods. Output can be obtained after each reduced mass-flow iteration and 
also after the final transonic solution. 

For all these variables except ISUPER, the integer given indicates the multiple of 
outer iterations at which the user wishes the output associated with the variables to be 
printed or plotted. A zero in any of these indicates that the output associated with that 
variable is to be omitted. A 1 will print or plot the output on every iteration. A 3, 
for example, would give the output on every third iteration. Any nonzero integer will 
cause output to be given on the first and last iterations and after the transonic solution, 
in addition to the output called for at other iterations. A large integer will obviously 
give only the first iteration and the final converged iteration. Care should be used not to 
call for more output than is really useful. The following list gives the output associated 
with each of these variables: 

IMESH Major output at every mesh point of the orthogonal mesh (i. e. , mesh point 
indices and coordinates, stream-function value, relative velocity compo- 
nents and total relative velocity, critical velocity ratio, and flow angles). 
This output is usually not requested, except for debugging purposes, 
since output called for by ISLINE or ISTATL is preferable in most cases. 

ISLINE Major output along streamlines (indicated by values in FLFR) at each point 
where streamlines are crossed by vertical mesh lines of the orthogonal 
mesh. Output includes z, r, and m streamline coordinates, relative 
velocity components and total relative velocity, critical velocity ratio, 
flow angles, and streamline curvature. Where the streamline passes 
within the blade region, an estimate of suction- and pressure-surface 
velocities is also printed. 

ISTATL Major output along station lines from hub to shroud (at locations specified by 
ZHST and ZTST arrays and at values of stream function specified by FLFR 
array). Output corresponds to that given for ISLINE with the addition that 
FLFR stream -function values are printed. 
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IPLOT 


Plotting indicator requesting output to be plotted on microfilm. Any nonzero 
value in IPLOT will cause the input data and generated orthogonal mesh to 
be plotted. Also, at each outer iteration which is a multiple of IPLOT, 
streamlines will be plotted, and meridional and surface velocities will be 
plotted for each streamline value. These will also be plotted after the 
final transonic solution. 

ISUPER Integer (0 or 1) indicating which solution (subsonic or supersonic) of the 

velocity-gradient equation is desired. If ISUPER = 0, only the subsonic 
solution will be printed. If ISUPER = 1, both subsonic and supersonic 
solutions will be printed. 

ITSON Integer indicating when information is desired for use in calculating input for 
the TSONIC program (ref. 5). If ITSON = 0, no information will be given 
for TSONIC. Otherwise, information will be listed with other output. 
Usually, the user only wants TSONIC data after subsonic convergence is 
reached. Using a high value for ITSON (i.e. , >20) will achieve this result. 

IDEBUG Integer indicating whether additional debug output is desired. If IDEBUG - 0, 
no extra output is printed. If IDEBUG > 0, the coefficients of the finite- 
difference equations, as well as 21 arrays of debug output on the orthogonal 
mesh are printed. Fourteen of these arrays change and are reprinted after 
each outer interation of the reduced- mass-flow solution which is a multiple 
of IDEBUG. 

Special Instructions for Preparing Input 

It is unusual to have no errors in input to MERIDL the first time any new data set is 
run. Therefore, input should be checked thoroughly before it is submitted. Errors are 
commonly made for the following reasons: inconsistent units; improper sign on u>, Vq, 
or whirl; input for arrays not agreeing with the input bounds for those arrays; and up- 
stream and downstream input not being of the form specified by LSFR, LTPL, and 
LAMVT. Also geometrical input into the hub and shroud arrays and the blade geometry 
arrays should be smooth enough that the hub, shroud, and blade sections will be fit well 
with cubic spline curves (see section (g), p. 22). These input geometrical arrays are 
all plotted on microfilm by the program after spline fitting is completed; the microfilm 
output will indicate whether or not the input was smooth. All output should be checked, 
especially from a new input data set, to see if it is reasonable. 

(a) Units of measurement . - The International System of Units (ref. 14) is used 
throughout this report. However, the program does not use'any constants which depend 
on the system of units being used. Therefore, any consistent set of units may be used 
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in preparing input for the program. For example, if force, length, temperature, and 
time are chosen independently, mass units are obtained from Force = Mass x Accelera- 
tion. The gas constant R must then have the units of (Force x Length)/(Mass x Tem- 
perature). Density is mass per unit volume, and mass flow is mass per unit time. 

Output then gives velocity in the chosen units of length per unit time. Since any con- 
sistent set of units can be employed, the output is not labeled with any units. 

(b) Damping factors FNEW and DNEW . - The input variables FNEW and DNEW are 
used as damping factors on the F f and 3(rV 0 )/3r terms of equation (Al). During each 
outer iteration of the subsonic or finite-difference portion of the program, new values 
for F r and 3(rV 0 )/3r are calculated. The calculated changes to these terms are often 
so large that if the full change was accepted on each iteration the solution would diverge. 
Using values of FNEW and DNEW less than 1.0 allows only a portion of these changes 

to be used. The value of FNEW or DNEW is the fraction of the predicted change in 
F r or 3(rV^)/3r to be added to the previous values of these variables. Therefore, 
if F r or 3(rV^)/3r do not converge, FNEW or DNEW should be reduced. However, 
this will reduce the rate of convergence so that caution should be used that FNEW and 
DNEW are not made too small. Because of the reduced rate of convergence, VELTOL 
is automatically reduced to obtain the same accuracy as when FNEW and DNEW are 1. 0. 

Limited experience has shown that for subsonic flows, where the gradients of flow 
variables are not large, converged solutions can be obtained with FNEW and DNEW equal 
to 1. 0. For most practical flows, however, this is not the case and some damping is re- 
quired. Values of FNEW and DNEW in the range of 0. 5 are most commonly used, al- 
though at times values as low as 0. 2 have been used. (Note that FNEW and DNEW do not 
have to be equal to each other. ) 

The user will have to gain experience with FNEW and DNEW before he can use them 
effectively to maximize rates of convergence. Maximum, minimum, and change values 
for F r and 3(rVg)/3r (DVTHDR) are printed with the output on each iteration. The 
maximum relative change in velocity at any point is also printed; it should gradually ap- 
proach VELTOL for convergence. By observing these values, the user can learn the ef- 
fects of changes in FNEW and DNEW for different runs of the program. 

(c) Hub and shroud flow channel geometry . - The hub and shroud geometry are spec- 
ified in the ZHUB, RHUB and ZTIP, RTIP arrays. Both of these curves must have the 
same z-origin (usually the blade leading edge at the hub). These two arrays must extend 
far enough upstream of the blade leading edge and far enough downstream of the blade 
trailing edge to cover the upstream and downstream boundaries of the orthogonal mesh, 
as well as the upstream and downstream input stations where streamflow data are given. 
If they do not extend this far, they will be linearly extrapolated and an incorrect flow 
channel may result. Relatively few points are needed to describe these smooth surfaces 
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(5 to 10 is a typical range for NHUB and NTIP) in order to have the program calculate 
smooth, accurate spline fits of these surfaces (see fig. 6). 

If the user knows the amount of the boundary -layer blockage along the hub and shroud 
profiles, he should revise the hub and shroud arrays to include this blockage effect. In 
this situation the output will represent a more realistic flow condition. 

(d) Orthogonal mesh . - The number of orthogonal mesh lines is specified by MBI, 
MBO, MM, and MHT; and the positioning of the mesh is specified by ZOMIN, ZOMBI, 
ZOMBO, and ZOMOUT. These four z-coordinates, all located on the hub, must use the 
same z -origin as all other geometrical input. The ZOMBI and ZOMBO locations are 
usually close to the leading and trailing edges of the blade at the hub, although they do 
not have to correspond to these locations exactly (see fig. 8). Mesh- size spacing in the 
horizontal direction is established by the relation between the four z-coordinates (ZOMIN, 
etc. ) and the numbers of mesh spaces requested. For instance, MBO - MBI evenly 
spaced mesh spaces along the hub will be located in the distance given by ZOMBO - 
ZOMBI. Mesh spacing in the vertical direction is determined by locating MHT spaces 

in the hub-shroud distance. Usually, none of the vertical mesh lines follow the blade 
leading- or trailing-edge lines since the blade edges are not usually orthogonal to the 
hub or the shroud (see fig. 8). So some vertical mesh lines will pass in or out of the 
blade region as they progress from hub to shroud. A suggested number of mesh lines 
between ZOMBI and ZOMBO is 15 to 30 in order to cover the blade adequately, depending 
on blade geometry. 

(e) Upstream and downstream flow conditions. - Upstream flow conditions are given 
in the SFIN (or RADIN), TIP, PRIP, and LAMIN (or VTHIN) arrays, which are all of 
length NIN. Downstream conditions are given in the SFOUT (or RADOUT), PROP (or 
LOSOUT), and LAMOUT (or VTHOUT) arrays, which are all of length NOUT. Upstream 
and downstream flow conditions are used, along with the assumption of conservation of 
angular momentum along streamlines, to establish boundary values on the upstream and 
downstream boundaries of the orthogonal mesh. 

Ordinarily, the upstream and downstream flow conditions are given as a function of 
stream function (SFIN and SFOUT). In this case these input values apply at all points 
along streamlines upstream or downstream of the blade; and the ZHIN, ZTIN and 
ZHOUT, ZTOUT inputs are superfluous. Values of zero can be supplied for ZHIN and 
ZTIN and for ZHOUT and ZTOUT, or blank cards can be used in place of these variables. 
On the other hand, if input is given as a function of radius (RADIN and RADOUT), legit- 
imate values must be supplied for ZHIN, ZTIN and ZHOUT, ZTOUT. In this case the 
upstream conditions are given on a straight line which passes through the two points 
given by ZHIN on the hub and ZTIN on the shroud. Downstream conditions are given on 
a straight line which passes through ZHOUT on the hub and ZTOUT on the shroud. These 
lines may lie anywhere in the regions from the blade edges upstream and downstream to 
the boundaries of the orthogonal mesh (see figs. 7 and 8). 
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The arrays of upstream and downstream input do not necessarily have to extend all 
the way from hub to shroud or lie on radial lines. They will be linearly extrapolated to 
the hub and shroud if necessary by the program, should the user only give data in a por- 
tion of the flow channel. 

(f) Mean blade surface and thickness coordinates . - The blade shape is described 
from hub to shroud by the arrays ZBL, RBL, THBL, and TNBL, all of which are two- 
dimensional (see figs. 6 and 10). Each of these arrays has NBLPL blade sections or 
planes, with NPPP points in each of these planes. When giving data to each of these 
four arrays, start each new plane of data (NPPP points) at the beginning of a new card. 
All the ZBL data for all the planes are given, followed by all the RBL data, etc. The 
origin for z-coordinates for ZBL should be the same as that used for all other 

z -coordinate input arrays. 

The input 9 -coordinates in THBL are of the blade mean camber surface and not the 
mid-channel flow surface. The program, however, obtains a solution on the mid-channel 
flow surface. It obtains that flow surface by smoothly fairing the inlet and outlet flow 
angles onto the input blade surface at an appropriate distance back on the blades from the 
leading and trailing edges. This distance is a function of solidity (see appendix F). 

The input blade thicknesses in TNBL are normal to the blade-section mean camber 
line and lie on a surface of revolution cutting through the blade (see figs. 6 and 10). So, 
in general, these thicknesses lie on a curved line whose ends are at different radii and 
are almost, but not quite, normal to the blade surfaces (see fig. 10). Because these 
thicknesses are used only for calculating blockage, it makes little difference whether 
they are obtained from a conical blade section or from a cylindrical section. 

The first blade section given at the hub or the last one at the blade tip (in the ZBL, 
RBL, THBL, and TNBL arrays) does not necessarily have to conform to the hub or 
shroud profile. It can be given within the flow region, crossing the boundary, or com- 
pletely outside of the boundary (see fig. 6). Extrapolation or interpolation will be used 
when necessary to obtain blade data where the blades meet the hub and shroud profiles. 

The user should attempt to give smooth data for the z-, r-, 9-, and thickness ar- 
rays describing the blade. The program makes use of spline fits of these curves to com- 
pute gradients which are used throughout the program. The microfilm plots of the blade 
sections given as output will indicate whether input data for these arrays were accept- 
able. 

(g) How to specify points for spline curves . - All the input arrays are fit with cubic 
spline curves for the purposes of interpolating, calculating derivatives, integrating, or 
any other required calculation. A cubic spline curve is a piecewise cubic polynomial 
which expresses mathematically the shape taken by an idealized spline passing through 
the given points. Reference 15 describes the method used for determining the equation 
of the spline curve. Using this method, smooth curves can be specified accurately with 
a few points, usually not more than four or five. Curves with uneven places, dips, or 
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highly variable curvatures require more points and are more difficult to fit properly. 

As a guide, enough points should be specified so that a physical spline passing through 
these points would accurately follow the curve. The minimum number of points to follow 
the curve should be used, since the closer the spline points are, the greater the effect 
of an inaccuracy in a coordinate. 

(h) Requests for output data . - The seven variables IMESH, ISLINE, ISTATL, 

IP LOT, ISUPER, ITSON, and IDEBUG all request different portions of output from the 
program. The optional arrays ZHST, ZTST, and FLFR are also used in this regard. 
The user should be careful to request only the output needed. Usually, only one of the 
three major types of output (at the mesh points (IMESH), along streamlines (ISLINE), or 
along hub-shroud stations (ISTATL)) is needed, since ISLINE and ISTATL output are ob- 
tainecTfrom calculated mesh-point data by interpolation (see fig. 11 ). Most likely, 
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Figure 11. - Location of three major types of output. 


streamline output (ISLINE) or station-line output (ISTATL) would be chosen. The fre- 
quency at which output is requested is also significant; this is controlled by the values 
given to IMESH, ISLINE, and ISTATL. 

ZHST and ZTST are given only if ISTATL output is requested. In this case, output 
will be given along straight lines from hub to shroud which connect corresponding points 
in the ZHST and ZTST arrays. Through the values in these arrays, the user can con- 
trol exactly the locations at which his output is given. These lines do not have to be 
radial. Typically, output is requested (through ZHST and ZTST) at several upstream 
and downstream locations, at the blade leading and trailing edges, and at several sta- 
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tions within the blade at percentages of chord from the leading to the trailing edge. To 
obtain output at the leading and trailing edges, the values in ZHST and ZTST should 
correspond (within a tolerance of 1 percent of local chord) to the intersection points of 
the leading and trailing edges with the hub and shroud. 

FLFR need be given only if the user wishes to specify the values of streamline along 
which he wishes output (when ISLINE or ISTATL are used). If the user does not specify 
otherwise (using FLFR), the program will automatically give output at 10-percent 
streamlines. 

(i) Incompressible flow . - Incompressible -flow cases execute as well on MERIDL 
as compressible-flow cases. In fact, they should converge in two or three iterations. 

No special input is required for an incompressible case, except that REDFAC should be 
set to 1. 0 to avoid the transonic solution which would be redundant and less accurate. 

(j) Straight infinite cascade . - The program is primarily designed for circular, 
stationary, or rotating blade rows; but the input can be adopted to apply to a straight 
infinite cascade as well. Since the radius for such a cascade would be infinite, an arti- 
ficial convention must be adopted. The user should pick a large average radius for his 
cascade, in the neighborhood of r = 1000. Then, since the blade pitch P is known, 

an approximate number of blades to use can be calculated from NBL = (2nr)/P. This 
value must be rounded to an integer for use in the input. With this value of NBL and a 
large mean radius, all the r-coordinate arrays of input (RHUB, RTIP, RADIN, RADOUT, 
and RBL) can be established about the mean radius, and a value of MSFL can be calcu- 
lated. Using this artificial input at large radius, the program will obtain a solution 
using a pitch which varies very slightly from hub to shroud and, therefore, simulates 
almost exactly a straight infinite cascade. 

(k) Choosing a value for REDFAC . - If possible REDFAC should be 1.0. However, 
if there is locally supersonic flow, REDFAC must be less than 1. 0, usually between 

0. 5 and 0. 95. For the best accuracy the largest possible value of REDFAC should be 
used (see the section NUMERICAL EXAMPLE). When REDFAC is 1. 0, the full mass 
flow will be used in the finite -difference solution for the stream function, and no tran- 
sonic velocity-gradient calculation will be made. 


OUTPUT 

There are four different types of output generated by the MERIDL program: 

(1) Main output - controlled by the variables IMESH, ISLINE, ISTATL, and ISUPER 

(2) Debug output - controlled by IDEBUG 

(3) TSONIC information - controlled by ITSON 

(4) Plotted output - controlled by IPLOT 
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Most of this output is optional and is controlled by the final input card, as already de- 
scribed. 

The output requested by the variables IMESH, ISLINE, and ISTATL is essentially 
the same for any of the three variables but is given at different locations for the conven- 
ience of the user. IMESH output is given at the orthogonal mesh points along horizontal 
mesh lines, as indicated in figure 11. ISLINE output is given along streamlines where 
the streamlines are intersected by the vertical orthogonal mesh lines (fig. 11). ISTATL 
output is given from hub to shroud along station lines (fig. 11) where these lines are 
intersected by the streamlines. 

In the following sections, output is presented from the problem solved in the. section 
NUMERICAL EXAMPLE. Since the complete output would be lengthy, only the first few 
lines of each section of output are reproduced here. In many instances, output labels 
are simply internal variable names. 

The following three sections discuss the different sections of printed output, the 
plotted output, and all possible error messages. 


Printed Output 

Table I presents the printed output from the numerical example. Each section of 
this output has been numbered to correspond to the following description: 

(1) The first output is a listing of the input data. Variable names are used as labels, 
and the output corresponds to the input form (fig. 5). This output is listed for every run, 
regardless of the values given to IMESH, etc. 

(2) This output corresponds to IDEBUG. It has three principal sections, as the out- 
put indicates. The constant quantities are listed only once; while the coefficients of the 
matrix equation and the changing quantities are listed in between outer iterations, as 
called for by the value in IDEBUG. The output giving the changing quantities contains 
two different sets of variables. (All the IDEBUG variables are described more thor- 
oughly in part n of this report. ) 

(3) This output indicates the calculated value of overrelaxation factor, ORF, to be 
used in the successive overrelaxation solution of the finite-difference equations. This 
matrix equation for the stream function u is solved iteratively (this is referred to as 
the "inner iteration") during each of the "outer iterations" of the reduced-mass-flow 
solution. This output, ORF, is given automatically on each run of the program. 

(4) This output indicates the maximum value, the minimum value, and the maximum 

change in the values of DVTHDR (3(rV Q )/3r) and FR (F ) calculated for each outer itera- 

v r 

tion of the reduced-mass-flow solution. The two quantities (3(rVg)/3r) and F r are 
critical to the convergence of the reduced-mass-flow solution, and the input variables 
DNEW and FNEW are used to control the percentage of the calculated changes of these 


25 



SilMli *i*s - laim X31N1 - auiua aossdaanUJ ivixv - i 1 vi»j 


O 

O 

o 

c 

c 

o 

in 

• 

o 


O O' 

h- O' 

i n -r 

CO -• 

PM • 

h* CO 

rg m 

• CM 

o 


O' en 

* ^ 

pm .<n 

in m 

m O' 

M f 

O CD 

pn 
I 


o 

o 

g- 

a 

«o 

pg 

o 


o 


IM 

O' 

rg 


pn 

in 


cr 


\ 


o 
o 
o 
* o 
uj o 
z o 

u. in 

o 


o 

>4* 

CO 

g> 

CO 

m 

rg 

o 


O' 

O' 

g* 


® 

CO 

eg 


<M 

g> 

g- 

m 

r\l 

o 

m 


O sO 

r* • 

nj -< 

r- eg 

in rg 

in cr 

IM PM 

O 


O 

O' 

m 

o 

IM 

g- 


-j cr 
O O' 
h- cr 
-J cr ►_ 

Ul O' < 

> *• »- 
• in ^ 
O O 


O o o 

O eg O 

O' O' O 

m in m 

m >?• m 

— 1 g> 

pg -d pm 

• • * 

o o O 


o 

OJ 

o 

in 

« 

o 

pn 

o 


O' 

<E 

r- 

r- 

g- 

PM 

o 


O' 

O' 

g- 


CO 

cr 

fM 


O' 

• 

rg 

in 

en 


o 

«© 

CO 


m 

I 


o <0 m 

— • « in 

in —« >o 

g> rg pe 

rg rg • 

■j- cr <r 

rg im *r 

O I 


H 

tn 

a, 

i-i 

o 

3 

CU 

S 

c 


o 
o o 
< a- 

U. O' 

o 

Uj O' 
a o 

• 

o 


cc 

z 



O' O' 

O O' 

in g- 

o — 

cc t 

pg co 

eg co 

• rg 

o 


O' en 

• eg 

rg pn 

in en 

en p- 

O O' 

wt rg 

I 


O >0 h- 

r- • rg 

rg -* O 

o pm o 

o- pg • 

pg cr 'O 

pg pg g- 

O i 


C 

z 


2 O 

c rg 


co cr 
cc Z 3 gj 
JCM 
11 PJZ 
in oc 


3 g* 
o o 
x r- 
o >o 


c o 

i i 

UJ UJ 

o cr O 

oo-o 
C O' o 

O' nj Q. 

a o 

o o o 

>c •-* 

• • • 

o o o 


a 

IT. 

O 

m 

CD 

O 

PO 

• 

c 



c 

CD 

en 


rg 

o 


O' 

O' 

g- 


CL 

CD 

rg 


O' 


g- 

g- 

g- 

r- 

• 

m 

pg 

I 


o o cr cr 

in —« O' cr 

m n g' g' 

pn -d — 1 

mm • • 

O «-' cc OC 

— en ® tc 

• • eg rg 

o o 


O' O' g a 
•- • •— 
rg rg cc «o 

in m eg o 

m rr C p- 

% • 

o o -* in 

— * pm g 

I I 



pm rg g- 

* w* w* 

O I 


© O o 'O >0 p^ 

me • • r- a 

pn p*' f-< ^ cr — ' 

n- p" ig pg m sr 

o n* rg rg ■ • 

o o o - cr pn >c 

-h m pg rg m g - 

• • • _< -h 

uo i | 


n- > 
a: 3 — s • 
ni g< 


> > 

«i < 

Cfi-C-d 
I a j a 
U' < UJ < 


a O a 
cc | O' 

< U.' 


« 3 O *— a 
a < *-' a — < u • 

u- c g 3 — i/s 


o z e o o o pg 

2 •- c a- o it- o a o- a «»-• 

<3 im — oxcr^-or-gj 

c;r‘ 2 g’ir“r-g , K»»' 

k mrv — i a h s h a n 




CJ 


e 


< *- r 

r- k i, 

« PK. • 


0 

1 


pr a 

r- o 
”* rg 
• « 
c o 


O' o 
O' O' 

g 


CD CC 
CO OC 

rg cm 


pg rg O p«. 

mm pr o 

pn pn — »o 

r-< «-i • • 

CO m ro 


Z o a 


a 

o o < 


<X 2 
Uj «“ 

a x 

h K 


UJ 

O wo 

O 2 <n O 

o <o «o 

h c: opio 
n <i o gj 
maif'c r 
*- •- re 

• • * 

O 


O 0 O' O' C pg 

O' O' • • 2 O O 

g'g'C. pgrg*—^g^ 

— 1-1 IT. IT. 2 no 
• • a n pr. p-g^g- 
(c<r a NH>a • 
CC ct O o • 

rg pg — • —* m g 

i l 


<1 K < 
K w * 
<l C. P 


o o 
rg o 
o n 


pg pg 
pg rg 

a o' 

cv rg 


a 2 
c 

I/. 3 


O K O O 

o 2 m c 

O' c« r- c o- 

c c «o c 

r- *c g- c a 

m a o a ft 

-d — 1 c\ 

C C c 


<o gi (- i/i r> 
• * 2 > — o 

— — * C. r»- ol 

pv rg 2 r r - — 

iv N n • • 

O o > cc — 

<g pg g 

i • 


26 


Ft AOE *FAN L IN 6 AND THICKNESS INPUT DATA 

TPL ARPAV 

f.OOOOOOO 0. 1 128000E-02 C. 63 07997E-02 0. I 36O200E-01 0 .22 2^900 F-0 l 0 . 33270 OOE-O 1 0 . 4442 200E- 0 1 0. 5 570000E - 01 


CO _« 


t 

U) 

O' O O' 
p- O m 
o* co 
a c p- © 

rg m rg 
«T -* 

-4 «$■ — « 


< i < 

UJ UJ UJ 

o o o o o o 

C) — < O m O O' 
g (T - ^ >0 N 
aomoo-4 
P-rgorgn — 
.g- -- in in 
»* -« 


OOOOOOOOOOOO 

0'000©«j-0®000in 


inOOO-OooO®0 


OoO'OOOoOOOOO'O'OoOOOvOOOO 

o®^®orJoO®»*u , '0®oo«-OpgrgcD®^ 

OrtOa'0'Na3 u >fno"J , (>(n'ONO‘r ta f\jMO«'t 

O' u> O' m O © ro^coco^fMr-^ocDroomogrw^ 
in O' ^-HfviPsj'^<4'»ntr>h-'OcocoO'0'oo 


o o o o o 

t ( < C t 

LU UJ LU Ui uj 

•O O 00 O O' 

O' O O' o O' 

® pi in ©o o 

p. f\j O <M 

O' CM ® O in 

NNO<V^ 
® O ® eP P“ 


< < 

Uj U.I 

o P- o O' 
o o O' 

00<J-N 
O 'J' CO O' 

ro 

—* Pi ® 


uj 

O' r- o 

p- O' m 
O' o <f 
r- r\j cc 
P- 'T ® 
'f sO in 

o -- 


oOOOOoOOoOOqOOOOoOOoO oooooooooooooooooooooo OOOOOOOOOOOO 


I 


I 


UJ 


I 


I 


I 


I 


I 


0 

1 


OOOOOOOOO'CJOOCJOOOOOOOO 

mo--ucMot rt - J 'to«^o‘MO'004'Oino'0 

v j-r--<f.- 4 ac* 0 'OOcr]'Oo'vj-cDmr— —<«*r*-'ToO' 

o-n^OtL^Ncr mPNirONPP'j'J *^in(r 
^n^'f^^^in^lAO'OO'OONOf'OffiO 


OOO'OOOOOO'OOO'tJO'OOOO^’OOO' 

0 P-'*pi kr *Cfn»npimOinn'*' 0 * T >'J’ 0 , fP“ -, 's;o 

^ir-^M7'O^Hcr'Oo<ocDm- - ^^a?ffl'f'0 

m®<t-p>»rs*.-o®oinr».in®o*op“0'OcocD 

® *n ® rn >0 O - pi 'T i£Di r 'Mr'tno®fnO'^ f 'J®'J' 
in 'O s cc (T ff i-i^rgf'j'P^ininh-'OcDxcJ’C'OO 
— ^.-‘.-(rgrgpjrgcgf^rorgrgrgrgrgrgrgmPl 


OOOOO O O O 

I » I 

UJ LU U-’ LU UJ UJ UJ UJ 

tnooooooooooo 
vff'ff'OOecOpomo- • 
CCO'^OCD^W'finMNP 
o O' ^ ^ in rg p- p- o *y .g- rg 
'T'0^>rmg50r'0h-0' c f' 
•t m ^ mj! - - cnf' -f ciin 
N'Ocoin-iin-c>f>-^-- 


OOCOOOOOOOOOOOOOOOOOO OOOOOOOOOOCOOOOOOOOOOO OOOOOOOOOOOO 


0 0 . 0*0 o o o c ococ? 

1 I I I I I I t till 

Lu Uj Uj UJ U. U. U- u_, UJ Uj UJ U. 

O O (7 OOcO'OOOoOOcOOoOlAOO' 
oop>ooor'" 0 -* cj 00 p', coorgcootr 
inmmO'M(j*co'tn f. in >o -j h m o; a in in 

o<ro' 0 '°** tf ' ,n '*-* p\org®.--oc ,r ' - oo 
orgo^o^omooor-cp-oxcooo® 

<rMf-<CM — r\j— <rM—ccv— ‘fMO'fvO' 


000*00000000 OOOOCOOOCCO 
t»cc«i^^«-ccro(J'00>-M fc( (c^ - ®fMinH 
m»r»ror w ^-o s cNia'P"^P' cc — c fr i cc O‘OcofMU^»-« 

inrep-a'OPgC'inr-ooNp-orgp- — 
r'-mf^f\jOO'f\j'#'®e£'( i n<viP-«©»-»®'^.- • o m cc m 
in'On-a'cro'^^N^’t'nninK'Oicttcj'O'Oo 
— — rgpgrgrgrjrgrgrgrjrjrgpurgcjpim 


O O 
t I 
Uj U. 
O P- 
O O' 
® ,0 
in ry 
V -4- 
in rsj 
in p- 


i i i 


i 


I 


I 


00*00000000 
ovcMOrJUincN 
^njino^ocoinnjin 
h -o nP CM^jOfn^-in 
'tCfyiP'PN^^inm 

m m ^ -4 n* fo cn vp a> ^ 

•'POpi— < rg~-.rg.^— --- 


OOOOOOOCOCCOOCOCOCOOO OOOOOOOOCOOOOOOOOOOOOO OOOOOOOOOOOO 


0 O o 

1 I I 

uj uj uj 
IT. O O 
O' o o 
CC o re 
P m ^ 
o in < c 

® c ai 

O -* © 

• • 


0 O o 

1 I I 

Uj UJ LU 

O cc o 
o o' o 

rg O' O' 

•-c o in 

m rg in 

«' co 'O 


0 C O c O o o 

1 I I I I I I 

Uj u; U; U. u.' Uj U; 

m c a c a op- 
O' O O' o O' o o 
rgooOp-'O'j- 
*v «p p ri co ® n 
P >o O P O P 
P P- 'll cc O W li' 
o -< O — O '- o 
• • • • 


0 o c o 

1 I I I 

LU UJ UJ UJ 

c m c a. 

o cn o O' 

O' P O M 
O M (C O' 
cc — • p. o 

U^ - Cl 

pg O' rg O' 


o r- 

o O' 

>P — 
n 
o — 

rg rg 
rg a <\J O' 


UJ UJ 
C •£ 
O O' 
P- O' 
cc in 
Pi rg 


i m - 


OOOOOOOOOOOOOOOO^OOOC 


I 

UJ 

OOOOO OOOOOOO'OCOO'O OoP' OO o 

rg ®®^ 3 >T'^--®nj--®CD*- J ®^ - »-*r^'Opjoc«o® o 

inN^o'ino^pPgjinmif'co'Oi’iffiin'O'^^'t — 4 

p-OM'fPO’tinPO'PtO'-'TOri'fa.inp^ 
ii^pMncN^ca pegr-,©.— o«P— •oP'-cum 
m«p«uirir»--jMfMtf^inir\PocuajirU'oo 
i^^-^^pgrgrgrgrgrgrgrgrgrgrgfMrjrgmm 

• •••••••••*••••••••••• 

UOCOOOOOOOOOCOOOOOOOOO 


0 o 

1 I 

UJ UJ 
o P- 

o a 
o- 

— O 
rg «o 


UJ 


0 

1 


I 

UJ 


00000*000 
OooO'O'rgOpg 
Of r lO' rn O'rgOg5 

— crorgoop-o' 

iTNOttaiTPa. 
cio^-om" rn 

fO cc nj >-i -j <-• p cvj 


rg i 


I 

UJ 

o o 

O *£> 
o rg 
rg o 
pj o 

pi <r 




CC°OOOoOOO 


o o 

t 


oooooooo 

I I I I I I I I 

Ul-UJujUJUjIUUJU. 
•occcr ® OR mo 
0'0' u '0'0'0'0'0 

O^ff 0*f «fC|R 
Pininmcgopo 
O'mO'ino'CD® — 
®r'.®®®0'cc»- 


0 o o 

1 » i 

LU U- UJ 

O' o o 
O' o O' 
'O o «• 

O P- CC 
—•mm 
cc rg p- 
® m a. 


0 o o 

1 I I 
U ' UJ IL 
O O' O 

o a o 
ci in R 
O nj © 
(O if p 
ci O «f 
^ cc — 1 


O o o o o o o 

t I I I I I I 


uj uj uj uj uj 
«Cor~0^ 
O' O O' C O' 

® NO 'C o o 

rv _. •— o «o 

cc O' O O' fg 

in in in <C 

® — OC QC 


O c 

O O 
m ^ 
p* p 


i i i 


J 


l 


i 


<J>0 O C C O O' C'OOCOoOOOO'OcO'OO 
©ininPin< , iNP'OONOinP'f©R©c’ OPJ i*' 
'KO' 'ONco'"- , omg , (rMOa'<Doo^®P - inN 
op p c <op rvPR©Ror>Rci^Po c? ooo 
■Cr. '0 ,J incDM , p®©ciNcc'C'CR'^'Cp fr Rin 
H'^P®(rO p '-HN«M^^UMnp^)ffCDo ooo 

— HMH^-NMNMMfgNMlNJN'MrgrgNr, K1 


m u uj _ _ _ _ 

OOO’OP-O'cOOOOO 
00'V'0'0'00* - 00'0*T 
OMr-cf'TciiPm [ 'inj 
pomajppinRRNc;^ 
fnin'OmP'Orr'J - ' - — < vr ''C 
®hP0CI' - NNPC i | 0C1 
—i 0 >©.— iin---j-^-^.~‘ro — 


OOOOO^OOOOCJOC CJOOUJOOOO CJOUJOCOOOOOOOOCJOOOOCOOO OOO 0 UO 0 O 0 c °0 


— pg 

000 

1 > I 
UJ in U’ 
-CO© 
O' o O' 
co© 

o- LT< 't 
m rg rg 
® pc cc 
P fNi p 


rg — rg 

0 O o O 

1 I I I 
© UJ u. UJ 
O o R © 
O O' O' O' 
0 4-0® 
rr o •"'* cr 
® C O IT 

o o p- 

© p in p 


CJ o 

I I 

II UJ 

^ o 
o o 
o in 
sJ- >T 


in r 


rg ^ rg w 

c o o C o 
I I I I I 

U. UJ UJ u LU 

® cc c in o 

O' O' O O O 

C in c in p 
in m o g- P 
n -o a - m 
O- o rg <o o 
p- r- cr p- •— 


0 o 

1 ( 

UJ UJ 

o o 
o o 
O o 
rg c 

sO <C 

in — 


000 

1 > I 
U uj Uj 
in OKI 
ooo 
O o « 
® o 00 
c a in 
in rg 'T 

-- p~ 


o o 

( i 

uj U- 

o m 
o o 
o o 
m m 
cc w 
m «t 


CCOOOOO'OOOOCOOOOOOOOOO 

® o om« 0 ' 0 rg^>m^c<iin®rg^^> 0 '^o 0 in '0 
« ^ in — o PNO"f»c®ttRi , 1 o«fcrffinP<t® 

mpsrrjnj-mOPc©0^'0©£)P-nj ,T 'NPJ 
m. rg»o^-<^®rgfncc®Pifn®«£;»-Lr»j'rgr w 'To'' iJ 
in <o p © c R-,jNN'f^ininp l o®©(rooo 
— — rgrgfvifvjP.rgrjfvrgpjrgfg'rgirvjmP’- 


rg *-• -- — . -- 

OOO O o o o 

III I I I I 

UJ U- UJ uj UJ UJ LU 

0*^0000 C O O O O O 
o O OrtOif, C^OgOo 

O©^(»i- j cx'-oog’inrvj 
cgpxpgj-opmco© 
pi ^cin^HO'O-oo'p-i 
goooo-p«-Hp.gN 
Pi O *-‘--rg^rd- J Pi^r r '-^ 


OOOOOOCOCOOCOOOOOOOCO COOOOOOOOOOOOOOOOOCOOC OOOOOOOOOOOO 


— rg — rg 

COCO 

i i i T 

UJ u u u 
P- o o o 
O O o c 
<r o p' o 
ginr- 
cooir 
p-m®r- 
O — < >0 rg 


O a 
p' a 
'C o 


© C 
I I 
U- U* 

g p 
o o 
o o 
o rg 


0 o 

1 I 
U' Ui 

>c c. 
o o 
n o 

gf O 
r*- yC 


— rg ^ rg 

0 O C o o 

1 I I I I 

u UJ Uj Uj u 

O O' *C O' ® 
OOOOO 
p - er a. 0 ® 
o — in a o 

r £i © © C 

•o m in © in 

m ff <5 r o 


OOO 

I I I 

UJ U- UJ 

coo 
Cc o 
nj in in 
-f in ^ 

— P'. rg 

-mm 
-• o ^ 


Uj U 
O C 
O O 


0 < 

1 a 
UJ or 
O <i 


: a. tt 
i «© a 

• o 


ooo 
rv» ® o 
tc o a 

in — >o 

in K)P 


COO 

in c rg 

c 'O 'T 

in rg p 
c m p 
a o o 


ooo 

P © ■# 

o©«f 
o r- P- 
rg m cc 
— — rg 
ro M\J 


C^COOOOCJCUOCOOCJOOOCJCO OOCOOOWOC 


OOoOOoOO 0-0000 

XNffl'lPlpK'O^C'f-CO 

©Ro^-PNPigo ® —< p- «g- 
cc«Nnj© a RffinNa rg p 
ttmpiao.-cr^'rgp'^ao 
niggininrcioocoo 

rgrgrgrgr-rgrgrjfvjrgrgpipi 

• »••••••••••• 

OOOOOOOCJ OOOOO 


J. o 
I- O 


0 c 

1 I 

uj a. 
3 P- C 

) O o 

> O «J- 
cv. rg 
® rg 
O 

o — 


i 

LU UJ U 


I 


O o 

o © 


c rv 
— rg 
• • 

o c 


O O 
c o 
Pi o 

— p- 

P-. p 

c o 

— rg 


OOOOO 
-rop-op- 
©i rvj -< o o 
-CCOO-J- 
g m g Min 

o «©• C P- o 

n •- n - 


o o o o © o o 


27 


TABLE I, - Continued. SAMPLE OUTPUT 



UJ UJ 

00 O 

O' O' 

4 J nO 

ip p- 

O eg 

<j © 

CD CO 

o o 


0 

1 

UJ 

O' 

o 


4 

o 

4- 

CD 

o 


<N 1 

O 

l 

UJ 

o 

© 

o 

c 


© 

o 


rg 

o 

« 

UJ 

o 

o 

o 

sO 

rg 

eg 

© 

o 


pg 

o 

< 

ui 

03 

o 

o 

© 

4 
in 


o 


eg 

0 

1 

UJ 

p- 

o 

o 

m 

o 

gj 

g> 

• 

o 


pg 

o 

» 

UJ 

o 

© 

© 

in 

© 

rg 

o 

# 

o 


pg 

o 

» 

UJ 

© 

© 

© 

© 

in 

r~ 

in 

o 


pg 

0 

1 

UJ 

© 

© 

4- 

© 

s0 

in 

• 

o 


pg 

o 

t 

UJ 

© 

© 

© 

4- 

r- 

<5 

© 

o 


o 

o 

o 

o 

o 

o 


o 


I 


Ul 


in o 
© o 
© © 
© rg 
g? p* 
•4- >0 


00 © © O 
©--O'© 
o CO f- 00 
in oo — i © 
m ip fp ® 


l 

UJ 


I 


mMMh - 


o o o < 
o 4 o < 

'f ffl lA ( 
4 — i rP < 
o rg pg i 


® • 


I'fiH'O- 


UJ UJ UJ UJ UJ UJ UJ (xj UJ UJ UJ 

P“OC0O©OP“O^“O(DO©oOO©OP-OOO 

0 ' 0 © 0 © 0 © 00 ' 0 © 0 © 000 © 0©000 

© O © o © O © O © 0 © 0 ©oOo© 0©000 

an-'O^OPONOffiOfflo^aNooooo 

OoNONO-ONOHOcDONOHO^OnO 

» 00 'OAO(TO^OhO®OPO'OOi/iOiAO 


d 

0. 

0. 

o 

o 

o 

o 

o 

• 

o 

• 

o 

o 

o 

o d 

o 

* 

o 

o 

• 

o 

• • 

o o 

d 

d 

• 

o 

• 

o 

o 

t • 

o o 

• 

o 

o 

• 

o 

o 

t 

o 











PJ pg 

pg pg 


pg 

pg 

pg 

pg 

rg 

pg 

pj 

pg 


PJ pg 

pg 

pg 


rg 

pg 


o 


© 


o 


o 


o 


o o 

o 

o 

o 

o 

o 

? 

o o 

o 

o 

o 

o 

o 

o ° o 

o 

o 

o 

o 

UJ 


UJ 


UJ 


UJ 


UJ 


lil UJ 

UJ UJ 

Ul 

Ul 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

ui 

UJ 

UJ 

UJ 

UJ UJ 

UJ 


UJ 

Ul 

© 

O © 

o 

o 

o 

o 


o 

o 

r~ o 

p- 

o 

p- 

© 

o o 

p- o 

gj 

© 

p~ 

o 

© 

© 

© 

o 

CD 

o 

00 

o 


pg 


pg 

o 

© 

o 


o 

o 

© o 

© o 

© 

o 

o o 





© 

o 

© 

o © 

o 

© 

o 

© 

o 

4 

g> O 

ro 


m 

in 


4 

4 

© o 

© 

o 


o 

o o 

© 

o 

© 

o 

© 

o 

© 

o © 



o 


o 


in 

in 

O 

4- 

O 

r~ 

_< 

© 

00 

gj 

cc 

o 

fp 


pg 

oc O 

m 

pg 

o 

p- 


pg 

© 


gj © 

© 

© 


g? 

r~ 

g> 


pn 


CD 


.-4 

00 

pg 

o 

m 

gj 

pg 

cc 

m 

rg g> 


o 



© 

rg 


© JT 

p~ 

00 

g) 

r- 

g? 

•—< 

g> 

© 

r- 


r- 



p- 

CO 

o 

© 

pg p- 


g? 

pg 

© 

00 

© 

gj 

<r 

CP 

•a* 

pg 

pn 

— * rO 

o 

m 


pp 

NT 


4 

<-i 

m 


m 


pp 


cc 


r- 


g> 


g3 

*■* 

© 


© 


© 


© 

i-4 

© 


© 

** 

© 

•—* 

O 

© 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d 

o 

o 

d 

o 

O 

d 

o 

o 

© 

o 

o 

o 

o 

o 



. 








pg 


rg 

rg 


rg 

rg 

pg 

pg 

rg 

pg 

pg 

pg 

pg 

rg 


rg 

PJ 


pg 

pg 

pg 

© 


o 


o 


© 


© 


o 

o 

o 

o 

o 

o 

o 

a 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

U 


u 


UJ 


UJ 




UJ 

u 

UJ 

UJ 

UJ 

UJ 

u 

UJ 

u 

UJ 

u 

UJ 

UJ 

u_ 

UJ 

UJ 

UJ 

u 


UJ 

u 

in 

© 

o 

o 


o 

o 


o 

o 

b 

a 

o 

au 

© 

© 

o 

© 

o 

o 

o 

00 

o 

g> 

o 

© 

o 

© 

o 

03 

o 


o 

o 

in 

© 

m 

© 

4- 

o 

sO 

o 


© 

o 

© 

o 

© 

© 

© 

o 

o 

© 

© 

o 

© 

o 

© 

o 

O 

o 

© 

o 


o 

o 

<r 

pg 

o 

pp 

pj 


r- 

4 


© 


© 

o 

© 

o 


a 

o 

o 

© 

o 

© 

o 

© 

© 

o 

o 

© 

o 

© 

o 

f\l 

IT 


© 

gj 

—4 


o 

r\ 

© 

g> 

4 

m 

g> 

p- 

in 



g> 

rg 

m-t 

IP 

pg 

a. 


m 

—4 


© 

SO 


o 

r - 1 

© 

4- 

in 

pg 

o 





o 

© 

© 

p~ 

*4 

O 

N* 


© 

p- 

m 

p- 

g> 

© 

Jf 

gf 

g> 

o 

r* 

CO 

© 

® 

r- 

o 

in 


4 

r- 


p- 

m 


p- 

pg 

o 

© 

gj 

P- 


g> 

© 



n 

© 

pg 

sr 

pg 

m 

pg 

PP 

w* 


—> 





-* 

-* 


-g 



-o 

© 

gj 

pg 

in 

pg 

©1 

pg 


pg 

g 

pj 

•g- 

pg 

•a- 

rg 

^a- 

pg 

<a- 

pg 

>a- 

pg 

o 

d 

o 

o 

o 

o 

o 

o 

© 

d 

o 

d 

© 

© 

o 

o 

o 

d 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

d 

o 

o 

d 


pg 


PJ 


pg 







pg 

pj 

pg 

pg 

pg 

rg 


pg 

rg 


pg 

pg 

Pg 

pg 

pg 

PJ 

pg 


eg 

pg 


PJ 

o 


o 


O 


© 


o 



o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

? 

© 

o 

© 

o 

Ul 


UJ 




Ui 


UJ 



UI 

UJ 

UJ 

Ui 

UJ 

UJ 

UJ 

u- 

UJ 

Ul 

Ul 

Ul 

Ui 

UJ 

u 

UJ 

u 

Ul 

UJ 

UJ 

UJ 

UJ 


© 

© 

o 

CL 

© 

© 

O 


O 


© 

© 

c 

o- 

p* 

© 

© 


© 


CJ 

o 


c 

© 

© 

© 

o 

o 

o 

© 

o 


© 

© 


© 

© 

© PJ 


© 



© 

o 


© 

© 

© 


© 


o 

o 

O 

o 

© 

o 

o 

o 

© 

o 

© 

© 

c 

-«T 

© 


a 

© 

© 

pg 


gj 


© 

© 

o 

© 

© 

© 

© 

© 

© 


o 

o 

© 

o 

o 

o 

o 

o 

o 

o 

© 

o 


PI 

m 


cv 

«a- 

o 

m 


— j 


© 

■4- 

pg 

IP 

fp 

in 

r- 


*a- 

g> 

gp 

PJ 

m 

© 

CP 

© 

pg 

o 

© 

o 

P- 

gj 



m 


cu 

© 


o 




o 


pp 

o 

gj 

4- 

p- 

r- 

in 

in 

© 

0 

in 

gj 

gj 


o 

eg 

gj 

a 

© 

r- 


© 

rg 

© 


u\ 

*4- 

gj 


g> 


pp 

© 

© 

in 

'T 




© 

gj 

r- 

4 1 

g> 

pn 

in 

rg 

© 

PC 

< 

«-» 

4- 

—• 



© 

<-4 

gj 

*-< 

g> 

-g 




© 

* 

<r 

'J- 

4 

•4- 


pp 

m 

pp 

pp 

pp 

CP 

pp 

pp 

ip 

pp 

pn, 

<p 

pp 

PP 

pn 

• 

* 

• 

• 

• 

« 



• 

• 


• 

t 

• 

• 

• 

• 

• 

• 

• 

# 

• 

• 

• 

• 

• 

• 

• 

• 

t 

• 

• 

• 

© 

I 

o 

0 

1 

o 

0 

1 

© 

© 

i 

O 

© 

• 

o 


o 

O 

o 

o 

O 

O 

o 


© 

o 

o 

o 

o 

o 

o 

O 

o 

© 

o 

© 

O 

o 













PJ 

pg 


pg 

pg 

rg 

pm 

pg 

pg 

rg 


rg 

pg 

pg 

PJ 

rg 


eg 

rg 

pg 

pg 

o 


© 


o 


© 


c 



o 

o 

o 

o 

c 

? 

o 

o 

© 

o 

O 

o 

o 

o 

© 

o 

O 

O 

o 

<? 

o 

o 

UJ 


UJ 






UJ 



U' 

UJ 

UJ 

u 

Ui 

U’ 

UJ 

iii 

UJ 

Ul 

iL 

UJ 

UJ 

u 

U‘ 

UJ 

u* 

UJ 

UJ 

UJ 


UJ 

o 

o 


© 

o 


o 

o 


© 



r~ 

o 

© 

© 

P- 



© 


O 

o 


p- 


r— 

© 

r- 

o 

P- 

© 

f- 

o 

«— 4 

o 

p- 

a 


© 

© 


pg 


o 

© 

o 

© 

© 

© 

© 

© 

o 

© 

o 

© 

o 

© 

© 

© 

o 

© 

o 

© 

© 

© 

© 

fp 

pp 

pp 

cc 

>4- 

1-4 

gp 

pp 

gj 


o 

© 


© 

o 


o 

© 

© 

© 

o 

© 

© 

© 

o 

© 

© 

© 

o 

© 

o 

© 

© 

r- 


>4- 

pg 

■4 1 

P- 

© 


© 


© 

—4 

o 

in 

-a- 

•a- 

g) 

r- 

in 

•4- 

© 

© 

pn 

*-< 

pg 

© 


© 

© 

pp 

pg 

© 

o 

© 



o 


© 

© 


U' 


pg 


© 

PJ 

pp 

pp 

pp 

in 

© 


g> 

pp 

© 


PC 



o 

in 


© 

© 

Pv 

pp 

pp 


*a- 



<r 

pn 

gr 


© 

© 


V 

pj 

't 

© 

o 

© 



in 

•O 

pp 

gj 

pg 

in 


in 

—4 

© 

1-4 

pg 

— < 

pg 


pg 


pg 


pg 

^-4 


pp 

g> 

pp 

in 

pp 

in 


in 

rg 

■4- 

pg 


pg 

4- 

rg 

•4* 

pg 

4 

pj 


rg 

't 

• 

» 


♦ 




» 


• 




• 




• 


• 

• 

• 

• 

• 

• 

• 


• 

• 

• 



t 

© 

o 

© 

1 

© 

0 

1 

© 

© 

1 

o 

© 

\ 

© 


© 

© 

o 

© 

© 

o 

o 

o 

o 

© 

o 

© 

o 

O 

© 

O 

o 

© 

o 

© 

© 

o 

—4 











pg 

PC 


pg 

pg 


eg 

pg 

rg 

PJ 

pg 

pg 

pg 

PJ 

pg 

PJ 


pg 

eg 


eg 

pg 

o 


O 


o 


© 


e 



© 

© 

o 

© 

o 

o 

o 

© 

O 

O 

o 

© 

o 

o 

u 

o 

© 

© 

o 

O 

o 

o 



UJ 


UJ 







Ui 


UJ 

Ui 

u. 

UJ 

UJ 

UJ 

U‘ 

u 

U' 

UJ 

u 

Ul 

UJ 

U) 

Ui 


UJ 

UJ 


U’ 

c 

o 

© 

o 

o 

© 

O 


o 

o 


© 

© 

© 

© 

c 

© 

© 

p- 

© 

p- 

o 

© 

o 


© 


o 

P* 



o 

o 

© 

pg 

© 

sa- 

o 


O 


o 

o 


o 

© 

o 

© 

o 

© 

c 

© 

o 

© 

o 

© 

o 


o 

© 

o 

© 


© 

o 


rg 

© 


© 

pg 

pp 

H 


pp 

rvj 


o 

© 

o 

© 

o 

© 

o 

© 

o 

© 

o 

© 

o 

© 

o 

© 

o 

© 


© 

o 


o 

© 

_l 

a 



© 


© 

■“4 


^a- 


IT 

o 

<r 


pp 

*0 

© 

fn 


** 

p- 

in 

o 

in 

r- 

4- 


© 


eg 

o 

pp 

p- 

<«T 


© 


o 

© 

© 


PVI 


© 

rv 

rg 

—4 

o 


o 

© 

pp 

© 

© 

© 

Ul 

<r 

CO 

•4 


© 



g? 

pg 


rg 

CO 

pg 

P“ 

pg 

g> 

pg 


gj 

o 

pp 

rg 

pg 

gj 

w 

*4 

c 

r- 

© 

in 


pg 

© 


© 

© 


© 

© 

© 

pp 


pp 


ro 


P' 

-- 

pp 



rg 

a 

pg 

P 

pg> 

o 

p. 

g> 

PJ 

© 


in 


in 


m 


in 


4- 

-4 

•t 

c 

i 

c 

0 

1 

d 

C 

1 

© 

© 

o 

0 

1 

c 


© 

o 

o 

d 

O 

o 

© 

o 

o 

o 

o 

d o 

o 

o 

o 

c 

© 

© 

O 

© 

© 











>- 


pj 


PS. 


pg 


rg 


pg 


pg 




rg 


eg 


pg 



c 


c 


o 


c 


© 


ee. 


c 


O 


© 


c 


o 


© 


o 


c 


© 


o 


O 



u 


UJ 


u- 




<r 


u 


UJ 


u- 




u 


UI 


UJ 


i 


u. 


UJ 


Ui 

© 


o 

© 

c 

© 

o 


© 

c 


o 


o 

a 

O 

a 

© 

© 


O; 

o 

a 

c 

© 


© 

o 

© 

c 

© 

o 


c 


o 

o 

© 

c 

o 

o 

© 

e 


© 

o 

o 

© 

O 

© 

o 

V 


© 

o 

© 

© 

© 

c 

© 

o 

© 

o 

© 

© 

© 

o 


o 

o 


o 

g> 


—4 


© 

© 

o 

© 

0 

© 

© 

© 

0 

c 


© 

© 

o 

© 

o 

© 

o 

o 

o 

© 

© 

© 



c 

t~ 

rg 

p- 

<4- 


gj 

fp 

u 

c 

P: 

© 

© 

c 


© 

© 



o 

g? 

o 

•4- 

© 

</ 

o 

gj 

c 

eg 

© 


o 

© 

© 

U' 

Pi 

U '■ 

u 

©. 

© 

© 


© 

g: 

© 

© 

c 


© 

© 


© 

© 


o 

m 

© 


© 

r~ 

o 


© 

© 

u 


o 

c: 


o 

© 

o 

© 

© 

p- 

c- 


o 

© 

c 

pp 

o 

r- 

c 

fp 

© 

© 

o 

© 

o 


o 

in 

o 

© 

c 


pp. 


sj- 


^a- 


r- 

-* 

P: 



© 

CC 

© 

p- 

c 


c 

-0 

c 

g> 

o 

gj 

o 

IT, 

o 

© 

© 

in 

© 

© 

o 

© 

d 

c 

© 

© 

© 

d 

© 

© 

© 

© 


© 

c 

© 

c 

© 

d 

© c 

© 

© 

© 

© 

c 

© 

© 

c 

© 

o 

c 


© 

© 


I I I I I 


o 

o 

o 

o 

o 

o 

gj 

o 


o 

© 

o 

o 

o 

o 

in 

o 


o 

c 

c 

o 

o 

o 

sf 

o 


o 

© 

© 


o 


o 

o 

© 


tp 


o 


© 

o 

o 

o 

o 


o 

© 


CO 

pg 


o 


o 

in 

© 

in 

m 

o 

—4 

t 

© 


<L 


o 


c 

o 

o 

o 

o 

© 

PJ 


o 

o 

o 

o 

o 

o 


© 



u 

< 

QC 

u 


UJ 

z 


© 

I 

UJ 

<o 

© 

© 

© 

© 

© 

© 

o 


o 

o 

o 

o 

a 

o 

© 

« 

o 


C; > > 

■“ < «■* < 

»- • a: o a c 
< a l a | 
h<Ui<JUJ 

to o o 

© C 

he* - ® 

h V! tf) 

1 M“ fv 
t KIAKI^ 


IL a 
a q: 
K <1 


o 

o 

o 

c 

o 

o 

cc 


© c © 
t I 


V 


28 


nilTOIJT PRINT CONTROL DATA 

ISl IMF ISTATL I PLOT ISUPER ITS 0 N IDEBIJG 



orsi^Osorur-m 

mr»r*iOC' 0, CC 

»-» tr ^3 r- f- o' r; c 
jr-tnrvjOotT'-'J- 
a >u ^ r> o h « 
t/i ct n r o» <\j rvj cm <\j 

OoOOC-COO 


— 4 >0 o a o'OO'co 

a. oi rr, >.*• ir* ^ r>- r*- 

OO'O^CT'OO'C'O'CT' 

ccJdocooc! 


OOOOOOOO 

«/5000oOOO u 

1/500000000 

ccooooooo 

_j OOOOOOOO 
D OOCOO.OOO 

x •••••••• 

ty~> oooooooo 

LLi 


unn u> u' tfMmr 

0 o o O o O o 

1 I I I I » I 

UJ UJ LLJ LL' UJ LU LU 
fSJ J/3 r- n't 03 f\J o 

O CMC o (T l^l 

-T O' ^ S >0 
O W+ O' <C sQ Q- <V 

<© in in in m if* 
• •••••• 

o o o o o o o 



0 OOOOOOOO 
CKOOOOOOOO 
XCOOOOOOOO 
>- -1 OOOOOOOO 
CtH-OOCOOOOO 
CCJOOOCOOOC 

* I • 

OOOOOOOO 

UJ 

1 


ff 1 ^ h* N Ifl N • - 

7 o ^ O M N 

(/) ••■•••• 

o o O o o C c 

ct 

a 




29 


ITFRATION ft, MAXIMUM RELATIVE CHANGE IN VELOCITY * 0.4765F-02 


TABLE I. - Continued. SAMPLE OUTPUT 


uutrtro-u-cj'cr 

— - 
c. K' iTi u ' u ' m u ' mm 

J ••••••• • 

C Q 1 cr d- ^ u t t-‘ cu 

I i l I l I l l 


oooooooo 
oooooooo 
a oooooooo 
a oooooooo 
<100000000 
tnoooooooo 
• ••••••• 

OOOOOOOO 


uu — crt\jfno<r , ^o>r 

I _j^-tn(\ja.irif\jOcra 

«/■- c i •*«•«••• 

uj 2 o 

a < ru rv i «x —i 


OOOOOOoO 

OOOOOOOoO 

XOOOOOOOO 

a:oO00OC>O o 

oooooooo 0 

UJOOOOOO 0 0 

Q • « • • J « • I 

OOOOOOOO 


OOOOOOOO 

OOOOOOOO 

Q-OOOOOOOo 

XOOOOOOOO 

<00000000 


j uj - m o o (r o to w 
O' 

tC>> I ( I » » • • • 

a: 


OUJ47Nnio<f<O^W 

c o •••••••• 

uj aif. « ui w m 


ccino^h't^N 

>or--r^unnf\j^ 

<T m r* fV rf) IP m ^ 
I OOOOOOOO 


cl oomomr-,$->c<M 

0 or^^f^ininfsir- 

1 C inh-Mtc^mrvjH 

h- I & xnv\-#&&'3--T 
a a oooooooo 


mr-oocroo*-' 

h c tfi fv c n o n 

. a.' O <r r- ^ m -< r- 

NN^OOO'U'ff 

I • • I * • » * 

i i i < i I i i 


a •••••••• 

ur.iAjff ert^u>^- 

lut-crcc^r-ccat- 

^■vf'f^cr-*-'rsjf r i*rir 


< «*■ < 

0 o o 

1 i I 

Uj U-' UJ 
0(0 
r- in cu 

i n et, o 
CJ h- O' •*> 
*- o- <r — 1 

X N tvjlP 

• • • 

0 c o 

1 I I 


m lti v?- 

0 o <-> 

1 I I 
UJ UJ uj 
cc <r 

it ,}■ i n 
•O o ® 
o m n 
-* m o 

• • • 

0 o o 

1 


0 o 

1 I 
U. UJ 

o o 
co k- 
cc in 
r> r*- 

<VJ f*g 


c o 


«r ^>no 

0 O o o 

1 I I I 

UJ UJ UJ uj 

m m *n ** 
c^erc 
a cr — o' cr 

G OJ ^ fO M 

u- «o ri < m 
C *■ <\j IT. r- 

• • • • 

o o o o 


>r < ^ 

o o o o 

• i i i 

UJ Uu U- Uj 
— <-« cr <\j 
O >U f\j o 

s r -i r> <*-. 
>f « o 

cr n od — < 
n <n .$• nj 

• • • • 

0 o O o 

1 I 


ox# 

Cl to * 
X UJ « 

OX* 


> < « 
»- r * 
- o * 

o o # 
o o # 

-I X * 


K # 

•- Ut # 

o o * 

c # 

— I !/> # 

Uj K * 
> * * 
— # 

a o * 

£*: 

a I # 

UJ V) ♦ 
K 1 U* 


* « * # * 

* * * 

• * « 

# X # a # 

# o * « 

# _i * # 

# u # • « 

# m # o * 

# CO * 2 # 

# < ♦ * 

* a # ft 

* * 2 * 

* * O * 

* CJ * *- « 

* UJ ft »_ « 

* o * < * 

* Z? # QL # 

* C3 « UJ « 

* UJ • K « 

* tt # »~ * 

* * # 

♦ * * 


Z 

c 


-J — • 

ujcttf»ooocc»» , 0'fvjm 

>*-U-OttU'P'DO^-J 

ik 

< V • • ■ • • • I • 

«-*01»OOOOOOOO 


• • ^•<Mn*o>QDrvjint) 
_j_j«-‘0'r(^<\!rvittcrr- 

tu uj 2 ••••■••• 

o;> — o‘f^'Ocr- - < s Jfn< 

O— ^-J(N(MP«jrsj 
csifsjrjpgfurvfgrM 


O 

it fNj'Cin^'Or-cr'T 

<i ,—*K'iccr\.'inK~<rCDcc 

h j i ••••«•• • 

imhOKfft^m^h.cc 
— *< — ^**-« 


• •cK^oir-.KH 

O • •~mr- r .oir\f-vr- 
*-_jX • •••«•»• 

aiu jrK^Kic^oo 
cu > «— hacttooctro 1 


OOOOOOOO 
.00.000000 
JOOOOCOOC 
J00O00OOC 
:OOOCJOOOO 
OOOOOOOO 
• ••••«•• 
OOOOOOOO 


mcsivma ooir 

<c h. a>c r ffj' - -<o 

a inr-f^cinr-srr- 
x •##••••• 
i/- r h < ^3 a w c c 
xr^a.O'QCaccoU 


• ro ^ it- o k- <r 


0-0 o o o o © o 
©oocooco 

OOOCOOOO 

acooooooo 

UOOOOOOOO 
COOOOOOO 
• ••••••• 

OOOOOOOO 


coooococ 

tnococcooo 

COOOCCOCO 

XOOOOOOOO 

K-COOOOOOO 

COOCOOOOO 

• ••••••• 

OOOOOOOO 




cc 


X • oooooooo 

<0*^00000000 
UJZ0O00OOOO0 
OCX^OOCOOOOC 

u. •*••••»* 

OOOOOOOO 


cn 


< cj r-fM^«rfn^-crcD 
i«a<»occ^u’n'«K(c 
c. C ££HN ^ ir ^®^0 

CL o 

OOOOOOOO 


C0O00OCO 
11111(11 
• U.UU.UJUJU.UJUJ 

j u o a lO^cvOoz'O 

<a*rt^KCM>uaH 
— ino^-aiMr^ 
xc^m^^ciM^rcocrr. 

<jo < »t ^ ri m m p* 

• ••••••• 

OOOOOOOO 
I I I I I I I I 
h~ 

2 J -* 

C a — 

a 

1 L _ 

cl 

Uj L - 

a o 


ir\ 


30 


a 

3 OCOOCJOCJOOOO 

CO •— OOOOOOOOOOO 
• ja •(«<*••<!«• 
Wuj JPO OOOOOOOO o 
Ul > — ’ 


3 OOOOOOOO 

m • *•. 0000000 ° 

• j a. t •««••• • 

tnu>»ooo©ooo° 

W > _ 


ac 

3 00©00000000 

VO *-00000000000 
' J W t « 

I-UJIOOOOOOOOOOO 
U >~ 


Ul 

z vi# 
O til * 
Q. 2 # 
5 *- # 
O 3 * 

o * * 
<1 * 

> u> * 
H CL * 
►—ho 

O «n * 

o * 

JC* 
IU 2 « 

> u * 

3 «* 
#. <1 * 


* -X * 00 * 

* C, * * 

* -J * « 

* U. * * « 

* <n * o * 

* in * Z « 

* < * * 

* 5. * * 

* * Z « 

* * o * 

* a * — * 

* LD * H * 

* © * •‘i * 

#3*0.* 

* O * Ul # 

* Ui * h * 

« a * ►- # 


3 ooocoooO 

CO •••.©000000 0 

I J «/) < t I • t • t • 

KU4»00000000 

o > — 

3 

CO 


<>wiftvfcO>00^-.N 

ui o: ci O' H o o o n n 

a D N N 0"0 

ho**»*«*** # 
m — 

W I I | I I I I I 


jiu-mooo^n'or' 
u_j.«icc'Om>rin>or-cr 
• C H •••••••• 

JZ UJ- - --H - H 

uj < a 


Cuj<ONl*IO 1 fU3>0^ 

*-_ia-if»«x C cif' rs »ov > <r 

aca • *•••*•• 

Ui 2 j — OoDN'OlfiWI'J 

2 


-J — 

u/ca'OOca>fcr'vio 

>HC*eca'(rooHH 

• h 5-£?'0 > o«or'r«-r^r- 

K < V l * • * • * • • 

— C£ 3 OOOOOOOO 


'TMnO'OJNtf'O 
• • O'fMMNflOtTK 

J J « • * • I * • • • 

UlUi3<M^'00'*HrMW<' 

ro<Mf\jrgrs.rv(SifM 


<Z • •• n a; Mi'f^ um ® 
»-jl »•«••••• 
• iUKO«n«nri'crcc 
3>*— 


o • «— in. KHOinr-»rK 

* J X ••#•*#!• 

cx.ujur^^»rr^cccroo 
UJ > — H OC CC CD ff> CL- O' O' 

3 — — — 


C C LJ o o «-< 
I I I I I I 

» • U Ul UJ UJ UJ ID 

co o «T (Oin r 

•-a n r- <n 

a C 3 m «c *“■• o- «— •»■ 
LUC^r-fvjr-fvjcDn- 
30 «T ^ n (O m 

0 c © o o o 

1 I I I I I 


0 o 

1 I 
Ui Ui 

tn sr 
o o 
h o 

OC 'J- 


0 o 

1 I 


<3 C r-fvJ^-^J-ro^iO'CD 

»-a*»occ-jar | UJN<c 
cc-a — » f\j >r in f- a. cr o 

coccoooc 


> 3: # 
O * 
UJ cl. * 
U * 
< UJ # 

u o * 
a < * 
3 -I • 

in ce * 

* 

O UJ * 
Z I * 
« H * 

* 

- T * 

in O * 
UJ 3 * 

-J C * 
C a * 

z x * 

<1 H * 
# 

» CO * 
CO 2 * 

HO* 
2 * 
Ul H * 
Z < * 
C H * 
a co * 
3. # 

O VO * 
0 3* 

o * 
>«« 
Ha* 
- < ♦ 
o > * 
o * 
OK* 
Ul < * 


V) c * 

£ 2 : 

H in * 


• * 


* 3 * < 

* 5 * 

* -i * 

* u. « 

* VO # O * 

* m * z * 

* <i * * 

* 3 * * 

* * Z * 

* * c * 

* O * — * 

* ui * K * 

* o « <. * 

* 3 * a * 

* t_) « Ui * 

* Ui * K * 

* QC * •- * 


o c * 
OH# 
-J * 
uj a * 
> 3 * 

1 * 
a * 
c 3 * 
— o # 
a a * 
uj u « 

H * 
Z VO * 
•_ UJ # 

2 • 
• * 


— C * 
HZ* 

o c # 

z -» * 

r < * 


* • CO UJ 

<>>-C04'O'} , 'f'f'O(\jC0'0 

UltfQHNO"CMrtNin’ - ' J 
HO • * * . a) in c\ ro o\j in co 
— | | IOOOOOOO 

III! 

* 

o 

_juj — 

lLj<®^f1t>^NOoON 

•Oh* • *•*«•••• 

(u<aimnf r i<r<t^<r^inin 


I 


C;UJ<KK'0'00''J-— ‘rjCr^cr 
»— j j rg^r-^'fotrv^^O'curM 
a. o a - * • * 


. j n o mo m i 

l <J -« -4 


i o ^ 'T ir\ 


ujca'fa^^onMa ® ^ m 
« o-rT'i o — m 

# •HJ t ^W?^‘HKf fc ‘KHH0CCC 

K < 3 v - ••••••••••• 

w* >“aJDoOOOUoPOCP 

oc »- 

o 

* 

cj n a —< r- © k m o cc ■£ 

z •• ^OM>rna^a «}■ m a. 

z uJu:J<-oM' 0 ^'ONfflinNa 

U or^*-fvjfvjfvfvnri^Nrin«oii 


ir, z fn-iinH^CO^o^m 

<i •«Mn^C’M« , 'Hooocoino 

Cr H3X ••••••••••• 

~ iujhcor'f*® a; ®h^ai^(o 

Cj _i>a — (Min^noiwatro- 

OC UJ <— I — ^INJCNl 


a: • oOfMc fVJ C'K c, r s -Hcr 

3 cj •— ofnsr--HKr*. u , '0*-‘V' 

I *-J» * •••••*•••• 

auioo''J'- j aj>otfi'j'*r'r(o 


3 • OOoOOOOOoOO 

<u — 00000©00©00 
UiZ30000©000©©0 
oc 3 ^ o ^ n n < in >0 r a? O' o 
HU •••**•••••• 

m ouducquoc. uh 


0©0©©C°0©0 0 
I I I I I < I I 1 I I 

UlUlUiUJUiUiujUlUiUjUJ 

O^rvjmf^^omcr ^ m 

3 m*Trn<MC\-f\jnj<vif\j r ^<M 
wtrimtr»inu'^mif>iriir'if* 

^ •“* 9* ^ ^ ^ 


C o 
— oc 
a c 
Ui o 

1 o 


3 C 

< OC 
— C 
X C 

< CJ 


OCo c OCO©COC 

I I I I I < I I I I I 

OCOOOCCJOOOC 

I I t » I t I I I I I 

UiUiUlUiUiUiUiUJUJUJUi 

OOOOOOOOOOO 

►'NNA/MMlMMlMvaN 

in ir. mi m ir. in in ir. m u', ir 


C C O c c c c 

I I i I I I 


T?‘ 


0 o © o © © « 

1 I I I I t 

• u' u. u' m u uj i 
— iC C av 4 'm d c 

<ia — rc'S'K c fn o c 

*-Ckh — modi 

x c w tn — «o pj h not 

v u c ^ n n (v a. 

coocoo 


I I » » I 


© © 
I I 


<i_ <r <c it o in ^ D ^ c c a 

►- a. — in r- c n mv r. k szm — 

C C aCNN^tt'fff(''> cm 
< C w m r- o - (v <r in f- c c ^ 

a c w ^ ^ r\ rv n n n. M" 

©co©cooco©o 


VO 


31 


52.41 1.153 


TABLE I. - Continued. SAMPLE OUTPUT 


a 

co • — r~-~*r~>ccv~*r m \c f i<TrrsV~\ 

• ex. •••••••»••• 

IAUJ^I^'t)mO^O® a 'N cc 

ujj-OH-Mtyjri^-ru'oo 

a: 


3 ONhOQJBO'f ^W>0 

m *<— ^n^rgr-fv — *t ec in «© c*- 

♦ j a * * » , » ♦ » » » * • 

UJ>^C 5 CTC 7 ' cr, CDCCCCOOCC cu CO 

tx 


OOOOOOOOOOO 
I — OOWO^JOOOOCO 


~ iAvj-^pj^o«c^[ro(NJ 

i/"- .-*co f N'a-'”0'»'(7'fsi^ 

• _j«/ i. •»•«»•••• 

►■u.iiitr'ga-r-^i^or-r-r-fpicc 
O i\> rg ix <\i <\j i\j eg rg rg eg rv 


z> C'-‘h-'Ttr'r\ja'«N>cc'3'o 

l/~. »-tf\rg^,iu^4r»ut*.o , 'yf«- 

i-LL 2 i^- , inMrjf<ir-c r- r- 
U > ^ o cr cd a? cc ar- a' o (r (T a 

«M — ^ — ***--*^.»g — •-* 


CZ. 

I? COOCCJCCOOOO 

CO • — COOctJCOOOCJO 

• _!«/»• ••••••*••• 

MdiCUCoUOOUOUU 

g> :> _ 


• CO 

J a O ^ O' fNj^corO^iMvOLArvj 
'D'voaiCOoOiArtONO'^ 
- O • • • • • « >o n • rg in 

o *-« rg eg -* rg — « • • • • • 

>— I I I I iooiooo 

I I I 

c 

> 

• c>- ••••••••••• 

JZaj(\j'r'CC'— •m^'Or^O'C 


o 

• *- i 

r • co uj 

t>*-<>o- 4 rom^og!Nja'CrcDr^ 

i ^••••••••••« 

— OoOOoOOOOoO 

I 


“-•d 


— rorMh-h-^m^r^rg^m 

^ CONhO 1 

H * • * * I * • I • • t 

a: | i -g-HfvjrvjrgpOfn 


— o 

* H* | 

• «/) LU 

[>~tnM-fi^fo^r~cDOAJ^ 

• O •(Mf-fOifi'Uf'-COO-aOU'^ 

^OOOOOOOOOOO 


JUj-'-OlB^fMO-NOM/lNN 
. j<r«ja;^in-to^<ONNi'J 

• C^- •••••••*••• 

jzujO'NfneiOffiNiro'M^ 
j < a. i I r\j <N rj i’ 1 


c uj < ^•crof^'Or-W'U'ir. oo «o 

— ’ _j x ^■o^fomom-o^no' 

a c a •••••••• ••• 


x <r <i ^ 


i 


i I 


)iu<ifNJOQcinf^'Oh-ir>cro-* 

•JltOlNjH^^MEWP'O'in 

: C & •••*••••••• 

J^_J^OOOOOOOOOO 
<3 < f | I | I I I I J 


• p- 

Otu<tsO-r^r-cDf\jmvj-ro(»iN; 

« <f (fiom 

Q! «j; D. ••••••••••• 

UJZJOOOOOOOOOOO 
J < « I 


uj C a *o 'C ia — — 
> — cj^-.mfoinine^cT'O-.min 

*h-ih-r^r^r*-r-r-r^a'cc«toc 

*-(X3:OOOOOoOOOOo 


f^r*tn'tf r -r- g; r*- u cr 

. • focccrf'-* - O'Oa.r-f^o 

> - m r. r -r ^ 'f' lo o >u n 


XOa 'O^'Tl^O'O f*. nT <M O 

>«— om-uirccf'-r-h-coa o- o 

•— or ^oooooocoooo 
a *- 
o 

a r\i— , r- a ' rv Jf' r >'C>r>j-f , i 

• • in r*- r*- ^ id Oh- fvj f^iTi ^ 

••••••••• 

Luajlomr^^-rnpor^'OcT'OCT' 

o:>wotrcrco<uii'ooa | ro(r‘r 

rg — *-*^~*-*^-*-+r+-* 


tuC“'Occccr^rgrgin--r-r-o 

>*- 0 'roocr-r-r-r s *eccoa;cu 

.►~3»©'Ommmmmmintf'iri 

K<>vi*»«» • •••*• 

_ a 30oooooc;oooo 
at — 
o 

COfVjCL <0 E CC' n «■- fPi g - in 
• • iBaco« , 'n^ - -J r 'Nfn 
-J _j *- • • • • 

a>-U(MCcccDa’o a. a- a> on 


C 

2 


z or- ^ cc di cr mpginr^ir- 
< • ~ oc^^rgco^^cc^f^co 
**•»•»•#»•• 
* u, h- o o u oaa-x>o>r^o 

_j>3 r 'j r r,mu>ir\OP-aU‘ c 


• r-f^o'-<o--f^in'rr-in 

c • --fg'Omrg,— r- in r- -• c 

►- j * *••••#•••• 

aLU30*-*«r»'Tooch »' <r 
uj > *" cr o » a: ® a) cc r> r r- 


c 

2 (ft'Tam'J'OlfiMEfOtC 

<3 «*~co | J^P r, (Njf r ^'- j r~ r ^mvj-o 

• u.'K(*iU(ro^ior--<NH«; 
-j > 3 ki N«rin^oochOo 

uj » i i —• 


• or^r-o^r-coOvoorj 
o • «* a <r •? i u 1 ^ h m m 
*-• .J X ••••••••••• 

ottu 3 fr, nf-rgoL-inor-m^m 
uj>-co s cco:r k r*r-'C^ , uir. 
x rg <~I ^4 — m4 


z ^}inr-'^'-^rg<rcr rv JO'- - 
<i • •-r-tn.-iomrgooccr-^gj 
hul*»i*«#'*»** 
*u-»-'(\<0'Or*-{\j'Ocoo^Of^ 
j > 3 m | ^rg^iP'OcU'OU 
UJ ^ I — — 

a 


h-a-— *moor^ccfoor->r 

o •i-uj^'C^(^a.''Ucc cc — 

«— _j » ••••••••«•• 

o:u.'3o^u'c*O rv -'0‘ o rn cr rg 
uj>wOO'a'fflr-r->o^'OuMi' 
5 rg 


X • OCCOOCOOCOC 

CO— ‘OOOOCCOOC'OO 
UJ Z XJOOOOOOOOOOO 

a r*'OH(Mn,tiPoM'ff o 

K U • •••••••••• 

CO OCOOCOOOCJO^- 

mrgpgcgcvjfv cg^-^— 
OOOOOOOOOOO 
I I I I t I • I I I I 

• • LUlUUjUJUJujUJLLilljUJLL' 

co fg»-«'Tcr'rr-or^Psjmr r . 

•-a^^^'J'fra^ctr.acf^ 
a C 5 rggjmr-f-m^orgmfpi 
LLC^orgr-rgr-tV'CC^cvro 
2 u fv >f ip r- as ~~ — , r-i ~- 

OOOOOCOOOOC 

I 

CCOOOCOOOOO 
I 1 I I I I I I ( I I 

• LU IX UJ UJ UJ U-> UJ Uj Uj UJ Uj 

_l c r. IT. ''-C'CC<rMC^'0 

■a a . — rag o h d f- r- r- 

>— c k o -J 1 ^ ip r- ir^ g M n r. 

xc ■w'Crgr-pgr-fgco-^rgpri 

o t_> 

CCOCCOCOCOO 


• oocooeooooo 

lO^OOCOOcOOOOO 

jzz:oocooc°oooo 

tD*'O' 4 W(*i<rm«Mc0'0 

- u_ ••••••••••• 

r, OOCOOOOOO^- - 


C C- ^JOCOOCO n >o n o c 

— a—rgrgcrOoa * 0 f'Ja>'O't 

acr a - fr >< r 'fNfM^ , 0' a: ^u>'Trri 

UJC'-'-‘*~ , ^*-i^COOOOO 


OOOOOOOOOOO 


_j c r-crom^crorgo*-*^ 
<a — r-crr^ccccorgainfr, 
^.r-K-rgrgrg r g^a'a.f^m>T f ^ 
Xc^— ‘-'^-^•-OOCOOO 

<lO — r* •-> -4 

CcCCOOOOOOO 


X • OOOOOOOOOOO 

< 0 «« 000000000 o 0 

IU2UOOOOOOOOOCO 
cx:Z«-0^n< r n^mg^r-aco'0 
k u. ••••••■••♦• 

is, OO^JOOOOOOO-^ 


oo noTMOir-mireco' 
•— tt«"»«o^<Nif\if\j~*fg<v<*grgfN 
a c j Kf^r-f-r-r^r-r-r~r-h- 

OOOOOOOOOOO 


_j cz & <x <r <x u O' O' o O' o o 

<a i^- 1 

»—CLh>-f*^r^r*r-r-r*-r-r“r~i < 'r- 
x C ~ m fp f 1- ! m ^ ^ M f ' 

•« o — — 

OOCCCOCOCCC 


oraj<t.’<J<cooHtjn 

.•Mf.Offnr >— fg cr cr 

c c q r. it. ip mt g- a fp ^ ct o 

r it. u ^ m ^ if. r- a cr 

- rg <g fg rg n. f\j rg n 


« C - - 
0.0 — ^ • 


• • • • 


OOOOOOOOOOO 


-J * 

< u u'CMNNCLri^rgp jo 

O o a m rs: o ^ o r. h o in h 

<c»'Oacrr.\ < fiPMao 

a u ^ f- — i\ tv m i\ i\ <\ rg r, 

OCOOCOCOOOO 


<j c r- <rr-og;a tro r^r-r 

► a- or. 

C(_a>j-fgO'ro ri o <v in r- 
•a c — ^cccr^-pg^ir. r cc u O 
a u -* —• — rg cv rg fg rg rg rg pj 

OOOOOOOOOOO 


32 


*»* INPUT OAT A FOR TSO w tC PR OCR A* •** 

****** 9 »************** 9 * ******** ******* 


U O 

a o 
o o 


rvj -r 
un <r 

<o ir> 


?? 

-s 


B? 

r- un 


0 o 

fU «« 

01 O 


I T | 

UJ UJ Ui 

♦ «M <0 O 

h- p 4 o + fM 
o (« o irn 
on® -f- 
a> 0 s ® r* m 
n o 9 1 
•- m r- -* •* 
# • • • • 
© © o © o 


♦ ff>h -Ow 
00^ -<lfl 

co r\ + v m 

OBNWf 

mm®®® 

t • • • • 

00000 


NNNNN 

00000 

1 T 1 » 1 

UJ UJ UJ UJ uu 
m 4 a* •* + 

0 << 0 (^-j 

<r m o >0 *> 
r* k — ® r* 

IN M O CO ® 
rg eg eg — « 
• • • • • 
OOOOO 


® ITS 

UN « 
V CO 
O O' 


*■ rN 
<0 'O 
f» £• 
sf O 
* ~ 


O O' 
f*N h- 
• • 
o o 


0 o 

1 • 
UJ UJ 

® *0 


? o o 
I ) 

UJ UJ UJ 
N M*> O 
— I f— o CO -o 
O ® ® m 4 * 
® **■ rg O' 

*4 rg r» 

• • • • # 
p o o o o 


co -4 r- ® r- 
•* O O **■ <M 
O' rg -4 ® -4 
r*- m -> ® m 
tr r- rsj rn •» 
«r m ® ® ® 


• • 


• • 


OoOOO 


N IN N N IN 

OOOOO 
I I • < I 
UJ UJ UJ UJ uj 
O -* -* eg ® 
m * a* © eg 

o ® ■— fw 
l*l<l'l-« 
O' O W ff> N 
eg -4 o ® ® 
M N N h *J 
* • • • t 
OOOOO 


X rg - 

<u -4 

© o 


r- o 

03 CD 
•** f- 


O' <M 
o ru 
® rg 
® m 


nr-M 
nr> o « 

N»N 0 

NWM 

ffl h- . 


o o o o o 


O (O ® <0 4 
tr oj r- cd 
o ® eg -- ® 
o ® r- t" rg 
® ® ** ® ^ 
m o *© ® 


00000 


rg rg rg jg rg 
OOOOO 
• I • • I 
UI uj UJ uj uj 

rg ** o ® «m 
r- o w in w 
O r- ® ® «o 

-* h- un o r- 
m ** o O' w 

rgNNHM 

OOOOO 


00 00 00 

t I >1 I I 

UJ UJ UJ Uj UJ UJ 

r- »-4 cr n ir\ c 1 o r*- »4 v 
r» cr -- rg o o 

rg eg ® p-»mmo in ini m 

® «»■ o o ® o o mr^<» 

® >© ® p n n z o uj®® 

®r --4 o in p _i • ® r- •-< 

rJ -O «« ®^ 4 ®a -4 — * ® -4 

• •• » • • «/l N •#• 

000 000 000 


3 O rg 
n O' ® 
1 ® rg 
3 in ® 
5 rg O 
n -4 43 
j r- <n 
• • • 


UJ Ui l 
r- -r e 
O m c 
rg ® - 
r- ® , 
O sT L 
® rg . 
rg ** ■. 


ITl®UN®P- 

® m r» r>». o 
eg ® it rg 
® r- rg m rg 
^ m ® ® 
«?■ UN 4 J ® 43 


ooooo 


rg rg rg rg rg 
OOOOO 
I I • I • 
Uj UJ UJ uj UJ 
® 43 «T -4 O 
rg®«ruN^- 
vr »r ® cr rg 
>o»r®Oh- 

® ^4 O O' ® 
rg rg rg •-* »-< 

OOOOO 


a. ® u. 
►- m a. 

E rg o 
rg t- 


I I 

ILI UJ 

* eg m 

m ® f- 
® •-* ® 
0 43 ® 

<r«' o 

N-H 

^- 0-4 


rg ® ® 

® ® r- 
rg ^ vr 
IN N r- IN 

me 
o rg o k 
un — r* ® 


?? 

UJ Uj 

* rg m 
® ® r- 
® — w 

O O m 

•T UN O 

rg — 


UJ UJ UI 
un «r -4 

® O' ® 


in tr o 
UN ® UN 
“4 O UN 


OOO OOO 


OOO CJOO 


OOO 
I I • 

UJ UJ UJ 
UN ® *$■ rg eg 
un ® r- o 
O' r- o'® -f 

■j-cg^-rg^- 

® o *r ^ un 

eg rg m o ^ 
® 0 «rN *4 — 

• • • • • 

OOOOO 


«f gj Cl MU 
m r- -a- ® ® 

a r> w * ~ 

un >r o n •r 
«j- un o o gj 


ooooo 


rg rg <g rg eg 

00000 

1 I I I I 

UJ UJ UJ UJ UJ 

r-r-r-rg^- 
® cc ® rg 
® c ® CD -4 

rg rg O O CD 
tt-MU'J N 
n rg — o' tt 
rg eg eg ■-* — 1 

OoOOO 


«j- a ® 
1 -4 o '»■ 

■ Ct C - 


a. •- 

it.hK 
n g ui 


r- O ® 
ru * 4 - rg 


® — ® UI 


000 000 


J l 
uj tu 
® 'i' rg 
o r* v 

r~ o cc 
rg ^ rg 
o ■# < 
rg un o 

® UN -4 

000 


un ® 

— 00 

m rg ^ 

rg o o 
«r o o 

woo 


0 q o o 

1 T 1 1 

uj u: ut uj 
o o c 1 ' *4 a 

»4 ® O ® »4 

— o^or- 
O "4 rg o co 
O ® ® -4 O 

P o 51 r 
® »»• cr -4 

• • * • • 

ooooo 


V o rg cc O 
rg "T ® — ® 
^ C PM 4 
<■ rg rg -4 
^ ® O ® < 
^ un ® ® ® 


ooooo 


K 


eg rg rg rg rg 

0 p o o o 

1 T I l I 

Uj ui uj uj 
^ 4 IfN ^ 

gj c rn cc 

g> r- o' ® ® 

® in o ® o 

fg^rgor- 
< rg «-4 O' « 

N N fN 1 ^ H 

O c o O O 


® 


® . 


-• — u O 

o o c O' ►- e 

<i , « IN z u 2 *— 

fw K -T UJ O' * C 

iuj >ca a' 
rg a O' * 


c r- o -- 'T cr 

— rg cr -4 00 

® m — ® ® >»■ rv. 

C O' ® N * 4 p- U 

-- o - tr a 

o o o coo 


UJ UJ u 

® ® u 
cc O f 
c c 
— <g c 
® r' - 
O u- 


tn o r- 

O' ® 43 

a. o W 


'TO 43 4 : 43 


U 1 UJ U' UJ 

un o o r- ® 

O ® "4 jn o 

rg -4 rg O o 

-- -4 o <7- rg 

m in rg 0 w 

rg in ® rg ® 

^ ^ *r o- -4 


COO C 3 o o OOOOO 


® in o un o 
— ® ^ O- rg 

rg ® cc o in 

a.' ffi 4 J ir C 
rg rg o in >r 
*r ® in 43 4 ) 


rg rg rg rg rg 

0 c C o 0 

1 I I I I 

LL< LL' UJ UJ U' 

5 ^ C M tt 

rg in O ® «T 

IfllNtl'f'- 
rgoo -4 — 

a- ® «r az a. 

•t rg O ® 

rg rg rg — -4 

C c O c c 


u 

C ® O a 

o -4 o O a 

c — r- < o 

a o <3 c a c 

■«j 1- i_ <r o o •— 

c g uj ■ iu o a 

• c in a o a 

® • 


< ~ — 1 a. 

a O o a C 

a 1 1 <3 » 

c lLa u> 

O C — 

O ~4 p> — «f «• ® 

— o rg o a O rg ® 
a o o cr v. ® — ® rg 
morgou >r in p ►- 
a o r ,iN|- 0 -ca 
— - — o *“ — 


^4 ^4 Q. •J' — -4 rg -4 _ -a 

CCaCCO>oOOc a 
1 I < r 1 1 < I » » I a 

Ui ui U U' UJ a ui LL! U: UJ < 

O O — ■ >T^CaCrfN~J;iNC r- r«j ^ 

O'-nN'ri'iM^r'NOMr 

- - 43 rg — * o eg o U 


rg O eg' O u cr mm. «crg~crrgoa*~o rg^a. 

a c o c v.r! r- u caraf^®v. cncr-o 

i7, o n c 1 c eg o. «r — o cr ^*®a — o <g® 
a ummi-OP ®a r~ tgr.uN® 4 >a •* tr- m <o ® 
o *r cr o in -o o r»t cc ^ — — 


UN o 

eg r 

NT ® 


c < o 


o *r cr o in 43 >r c ® a. 


a 1 » t 1 1 1 Jc o «3 

<3 ur in u: ur ur u. < 

o- -c O m 1/. in <• 

v - 4 r-®® 4 >o u 

®q cr r- c o m m »/■ 

r. n. ^ o v 4 r ^ r u 

«r 11. in U' u' ir a r-'Q — 

4 j a mig — o a. a 

— rg rv n ^ - _i 


ooooo c o c o c c o cooocc 


O 'ftt'fO* 4 £ l r\l 

<C< 3 C 0 rC<c* 
I » < t l 

Ut U U. UJ 

c a m a “’'Cl 

O a rg Q rg Q r q <f ( 

Cu.au.Um. cue i 
o_CoOou»-r“i 

o^-a h-a-t-OK ® 1 
O q. ® o. n a in a - 1 

o — — «■ «, 

o c c o c 


V. 


y 


33 




34 



zj oo’-’ooooo 
10 •-^oo'-'ooooo 

» _i O. **•#•*•• 

iouj.*OoCJOOOOO 

tu > ~ 


* ft =7 # 
* -J * 

• c * 


C3 * 

O ft 
X ft 
K- ft 
UJ ft 
» ft 
ft 


C * 
-J * 
u. ft 
(O * 


§ OOOOoOOO 

to • — oo°oo°00 

I J I/) I I t < I I ' t 

KUU300000000 

s*- 


a - • to 

<t>~io4-comsO'4.-*ej 

uiq:oO'^ h o-oonn 

►- U t •«••«•* • 

go 

— • l I • l I • • 


w 

a- _j<rm!nmmmr-coo 

• CH- ••••••*• 

..jZtJJw^-^w--— -jeg 

UJ<ir*M^^^wnfrtn 


CLU<0' [N J( <1 0'J'^'0 > f 

a c & ••*»••*• 

Mi2-j^oa3r*»oirfON 
a <x<ojrg^H^^»4^4«-< 


O 

z 


ox* 
a- <o * 


> < * 
t- z * 


* « o * 

* -I * 

* c * 

o « co * 
c * * 

X * 3 * 
K ft O * 
UJ ft _| ft 
X ft U- * 

* to * 


Z to * 
C w * 

o z « 
a *- * 
c -> * 
o a * 
« * 

> UJ * 
k- a » 
— h- ft 
o co ft 
C * 
-i C * 

uj z * 

> O * 


> * k- ft 

* I- * 

> * u * 

CC ft O * 

* -J * 

* uj « 


• CMTOOc^^O 

D t-'000<tf'J\ffiC 
•- J* •*»*»«»* 

auJicrfN^voajcro'o 

aj>c.^oc® e t®o:®o' 


UOUUOUUO 
I I I 1 I I I I 

• • UjlAJUJUJUJlljUJUJ 

O o ocn^encn^en,* 

►.a»-<^«inMvOo 
(t c i 

ujG~r*>f\jr-cgooncD** 
X UJ 'J’-tf'mmpgeg.-H.-* 

»»♦•»#•# 
OoOOOOOO 
I I I I I I I I 


* -J * 

z <■ * 

o * 

*- k- ft 
K « ft 


> * k- • 

* * 

> * O ft 

ec # a # 

* -j * 

* uj * 
• * * > * 

* * 

• -» ft 

* -j * 
« < * 


uj « y ^ r o 'f ^ -r 

1 j ► y (\j a: it m O O' 
COCJ • ••••••• 

uj 2 o Ho<r.f- <c- m f" eg 

2 


3 

O 


w .d 


-j'OOfvccU'J'fn 

<rinmin t/> >n* , -en© 


OUKJtfNOO-fO'O'f 

M ji^^oir(\joao 

c* e: o •••••••# 

UJ 2 J ^ O ^ ^ ^ucrifvi 


nidi oc <\i_«C (Tmcc — 1 ro 

>~-ur w ccO'ooo- - ^ 

V • • • • • I • • 

*— a.30o<~>coooo 


m >r gj . 
U> -4 >4 « 


* 

* 

« 


Uj 

« 

CO 

ft 

h- 






►H c 

ft 

ft 

ft 

* 


UJ 

ft 

CO 

ft 









* 


* 


K 

* 

*0. 

ft 

C_> 

• 





o o 

ft 

ft 


• 


k- 

ft 

< 

ft 

Uj 








* 


« 


<3 

« 

X’ 

ft 

z 

«J 





O c 

ft 

ft 


ft 


<1 

ft 

X 

* 

z 

• 

• 


O 

gj 

r- 

n 

* 


» 

z 

X 

* 


ft 

z> 

UJ o 

a fsj -* o 

cc m 

CC 

ri 

-J X 

ft 

ft 


* 

z 

X 

ft 


ft 

*— 

— J 

-J 


CC 

g> 

u\ 

— * 

# 


« 

o 

k» 

«■ 

Cr' 

ft 

u. 

> — 

UNCDtr 

O' o 

o 

** 

UJ h- 

ft 

ft 


» 

c 

*■* 

ft 

C 

ft 

•J 

u. 

Uj 

3 

• 

• 

• 

• 

« 

3 

* 


X 

« 

z 

* 


• h- 

X o <5 -O 

<£> r- 

r- r- 

n- 

> a. 

ft 

ft 

3 

* 


X 

ft 

z 

ft 


CL 

> 


o 

m 

gj 

O' 

* 

O 

« 

k- 

Vj 

« 

*— 

ft 

X. 

k- < 

s * • • 

• » 

• • 

• 

o 

ft 

ft 

C ft 

*- 

CJ 

ft 

•— 

ft 

X 




•M 

•— i 

»-• 

— < 

ft 

— J 

« 


a 

» 

3 

ft 

<s 

*- a: 

W 

0 

0 

0 

0 

0 

o o 

o 


ft 

ft 

_J 

» 


ex 

ft 

3 

ft 

CO 




eg 

eg 

eg 

eg 

* 

lo- 

* 

Zj 

CL 

•» 

l_> 

ft 

UJ 

a 





CO UJ 

ft 

ft 

u. 

« 

-J 

Q 

ft 

C 

ft 

UJ 








ft 

co 

♦ 

c 

a 

* 

X 

ft 

a 

u» 





UJ I 

ft 

ft 

CO 

« 

G 

O. 

ft 

X 

ft 

X 








* 

CO 

* 

to 

< 

O’ 

o 

ft' 

►- 






— ►— 

ft 

ft 

CO 

• 

to 

«3 

ft 

O 

ft 









* 

«a 

ft 



ft 


ft 

VI 


o -o r- 

n m 

ft o 

o 

►— 

ft 

ft 

<r 

ft 



ft 


ft 

— 1 

C' 







* 

X 

* 


k- 

ft 

z 

ft 


• • 

CO >o IT! 

-* o 

ft ft 


•— u 

ft 

ft 

X 

ft 


K 

ft 

z 

ft 

•a 

2 



eg 

'J- 

in 

O 

* 


* 

o 

z 

ft 

< 

ft 


-J 

— . • • 

• • 

• • 

• 

o c 

ft 

ft 


« 

o 

Z 

ft 

< 

ft 

z 

< 

« 


m 

OC 

eg 

in 

* 


* 


UJ 

ft 

X 

ft 

1 

UJ UJ 3 D o o 

O' -< 

cm m 

ft 

o 

ft 

ft 


ft 


UJ 

ft 

X 

ft 

o 

h- 

-J 

I 

• 

• 

• 

• 

« 


« 

z 


ft 

k- 

’ft 

1 

OC > 

— ■ _ i ^ m 

<N» 

rg cm eg 

-J CO 

ft 

ft 


« 

z 

•• 

ft 

h- 

ft 

cr 

• 

UJ 

k- 

o 

•—< 

eo 

■4- 

* 

-J 

« 

c 

o 

ft 


ft 



<v evj eg 

cnj eg 

eg eg 

eg 

UJ H- 

ft 

ft 

-J ft 

O 

Ui 

ft 


ft 

a 

— J 

> 

3 

*-< 


H 


* 

-J 

« 

to 

< 

ft 

a 

ft 







> z 

ft 

ft 

-J 

ft 

CO 


ft 

oc 

ft 

X 

Uj 


w- 

•-j 




* 

iJ 

* 

z 

CL 

ft 

UJ 

ft 


« 






ft 

ft 


* 

z 

a: 

ft 

UJ 

ft 

k 

Q. 







ft 

u. 

« 

<• 

o 

ft 

-J 

ft 


c 





a c 

ft 

ft 

u. 

« 

<1 

e.- 

ft 

-J 

ft 

a 








* 


«■ 

a 


ft 

— J 

ft 


z 

<\j 'T in 

•o 

n- ce 

ft 

C O- 

ft 

ft 


ft 

a 


ft 

— 1 

ft 

u 










ft 



ft 

«* 

ft 

a 

<L .• 

rh a: eg 

«> r- 

OC CO 

au 


ft 

ft 


• 

K 

> 

ft 

<1 

ft 









* 


ft 


►- 

ft 

X 

ft 

UJ 

k- _j 

X • • • 

• • 

• • 

• 

at x 

ft 

ft 


ft 


K 

ft 

X 

ft 


• 



O' 

\ n 

o 

O 

* 


ft 


►— 

ft 

CO 

ft 

cc 

• UJ K o ^ <n 

ft- U' 

o r- <c 

UJ CO 

ft 

ft 


ft 


►“ 

* 

CO 

ft 

«<a 

U> 

• 


«0 

o 

gj 

OC 

ft 

« 

* 


u 

* 


ft 

X 

-» > 

3 -• — 

pM 



k- UJ 

ft 

ft 


« 


O 

ft 


ft 

k- 


-J 

X 

• 

• 

• 

# 



* 


c 

ft 

to 

ft 


UJ 


M •— 1 



z X 

ft 



ft 


(J 

ft 

to 

ft 

z 

OC 

UJ 

3 

a 

CN 

■4- 

gs 



* 


— i 

ft 

Uj 

ft 

z 

(X 






ft 



* 


— J 

« 

UJ 

ft 

o 

UJ > 


h~ 

CL 

CC 

cc 


rvj fr, sr 
eg eg eg eg 


o 

X 


X • OOOOOCOO 
<lO* , ‘ t) t)OOOO u O 
UJZDUUOOOCOO 

a ^—-oocoocoo 
* **-••-“ 


k- a 
to 


00©00©C0 


<1 c 

►— a'"* o ®-'4 -u 'ro,0r-cr 

OC o h«-hh^h^h 
UOOOUCUO 


<c k eg ■4 •* m — o' cc . 
•- a -« O ec •J’ O n>or*'a:> 
OCa--fvj< , ir\f>a;ci o 
<C~444.4444cn 
ty y~) wK^***+^**~*w+ 

ooccococ 


co«->©ooow 
I I I I I I I I 
• U< Uj UJ UJ U, UJ U-' IU 

jc CTtJvrfvjoct) 

<sa— .rvj^-P'opO'Oa.^ 

<i o *4 >r 

• ••••••• 

OOCOOOCf 


OCOCJOCOC 
I I I I I I I I 

• u.' U- UJ UJ U- Ui U. U.' 

_j C o cl 2) -J" fv c; a' vO 

■«a q_ r\. »r r- c pi € o. - 

c x r-^-inC4<£^r- 
p>j ec 

< c >T>j , r‘ © eg eg 

• ••••••• 

OOOOOOOO 

I I I I t I I I 


. X ~ 

ij"""MO'irif. tr* a* 
. c — 


35 



'INC T I n N f INTFPICP VFICCITTES. VFirCITY CC w PONFNTS * ANGLES. AND SURFACE VELOCITIES 
ALnwr, LINES PPO« Kill Tn SMPrim AT VAPincS STATIONS thrdiich the plaoe row 


— ooooooooooo 

OO * — 00000000000 

•J& • • • 


“ OOOOOOOOOOO 

uO *-00000000000 
• _J (/>••••••••••• 

►-UJXOOOOOOOOOOO 

o > — 


• VI UJ 

t>*-®lT'Qnf'JCO'fON<3^ 
jajQ-HMo-'ON^srin^-^o 
tr;v.r^co-*fMmo'r~irTfvj(M-^ 
- O • • • i CD U> (\J n N If) (C- • 

I — *•••*•• —i 

— I | IOOOOOOO 

4 • ' 1 1 

5 

JlU'*NNl<lOcC(M'Ti/>ifi'Oo. 
. _j < O' vTfCiO'— rsjO^jOfMfvj 

• Ch- »••*••*•••• 

JZ^-inaomifiNcco^N 


OUJ<a:No^^'f ,H rsity^o 

— — ^crrsicsj 

or C Q- •••••••••*• 

uj z - 

X «I <J — — | I I I I 


■ 

ujcqc(m 0 it. 

> — o — p - o — r^' 3 - <jf^cr— <rr 

►- <i ^ • • •• **•••• * 

^-or^ooooooooooo 


_i-j— ••••*•••••• 

ujUJ 3 nctMO--<i(\;®inwff 

or>'-f\'(\ir\jrv.frpnvr N j-ir'vO'f 


C 

Z rna mc^tro^ciniri 

< •-iriiDO'M^Oiniririno 
»— -j X •••••••*••• 

• ujH-®r^r^oc®cLrwotn^fo 

_j > -*~*cgrT«rm'Uf*-a;o'o — 

UJ ~ ^ — 


c • — O' < h h a- fr, c *- i r- 

H J s ••••••*•••• 

or U-' •» O C->T— CL>«o»r\*yN?->^'T 

— ccr-r-»^o>Go>c«ooo 


3 • COOCOOCCCOO 

"aO^C'OOOOOOCOOC 

UJ Z. Z2 o o O O 0.0 o o o o o 

K u. •••••••••*• 

V, COOOOC^C'-’OO — 

CO0CO0CO0O0 
< I I I I I I I I I I 

• • UJ U5 UI UJ Ui IL Ul Uj UJ 1L Uj 

co c r ec c ^ - a — r*- no 

*— a— — — <^-r^>T'£j 

ac s in^^fsjrvfNjrvjfMPvjrvrj 

U-l c «— U 1 IT' 1/1 If! \X 1 Ifl i^ ir IT 1 X 1 IH 

3 O <-• < — — _l _4 •— > 

OCoOCCCOoOC 

I I t I I 

OOC.OOCOOCOC 

I I • I I I • I I I I 

• lL‘ U- UJ UJ 11 !!•' LL' LL IL U. U.' 

— c- e c* c c o c c c c o e 

<rar — — — — ~ — — 

x c » in it. u in in it- ir.ir ir in ir 


c o o c c- o c i 

i i i • i i * 


c C- UI M -LJ r Ui r, in U' o 

~ a — u-> r anv-fufir-. — 

C c crcrvfv — a ^ c r* r- c* <*■> 
« c —■ w' a — 1 f' >r r- b c - 

12 l_ ^ - 4- IN l\. M\(\ Mil (T 

ooccoocjcooc 


c re. C. a 

o f- U c > 

O — <3 C 

z o *. — u c 

^ O ►- >: D C - 

C- -T Uj • U.’ C 3 

• c <M a • 3 

— m «- 


36 


o o 

• I 

UJ UJ 

-• o r- m 

— < O CM sj- 

^ O' in O' 

cc o m 

cj cr m xy 

r- in o — < 

fn co — < -4 

• • • • 

o o o o 


0 o 

1 I 

UJ UJ 


o O 
• I 

UJ UJ 

o o 
eg m 
eo cr 
~H o 
p- ni 

O o 
in r- 


o o o 

( » » 

UJ UJ UJ 

r- «y rg >o O' 
r- o .y rg 

o r> « in r 

o r. «j h ^ 
O' O' cu r- c<i 
(TOITN^ 

-• m r~ ^ ^ 
• • • • • 

O o o O o 


jvOinK- 

^ o> r- *o m 
© o ^ hu> 
® a ^ a w 

O id rg fe «y 

^ >o <o 


o o o o o 


Csj f\J r\J <M <N 

0 o o O o 

1 I I I I 

UJ UJ LU UJ UJ 

O <© eg o O' 
rg in o m 
in x© ** O' «y 
<NJ O **l *0 CM 
a> >o ct 'O in 
eg rg CJ O' O' 

• • • • • 

O c ° o o 


0 o o 

1 I I 

u. UJ UJ 

xy sc r- f- 

»-« cr «J- m 

mm >o •© 
m co p- r- 

O CO "O' o 

o 0 m eg 

re r- O' ^ 

• • • • 

O o o O 


•— 1 .-4 

O O 

• I 
uj uj 

^ £ 
— < O' 

m ^ 
m ® 
O' ® 
o O' 
m n- 

• • 
o o 


0 o 

1 I 

UJ UJ 
CO © 


~ eg 
© 
m co 
m 

'J- r- 


o o 


O O CJ 

( I I 

UJ LU LU 
rgp-mop- 
© m ® eg o 
p- o co © 
o co co in *4 
O r~ e~ rg o 
ao m eg m 

WfijNW^ 

III*! 

O o o o o 


co-*r*-cor*- 
•h o* o r- eg 

O eg -- m 

f-. m —* ® m 
on*egr»'T 
»y in © © © 

-*4 "H ^ * 

• • ♦ • • 

o © o © o 


m cj eg rg eg 

0 o o o o 

1 I I I l 

UJ LU U> UJ UJ 

rg p- O' © in 

a>o^ H 0 
o in m p- o 

O eg m >t O' 

o ^y o p- m 

rtrvi*-CM> 
nj nj eg »*- — * 

• • • • • 

O O o o o 


0 o • o 

1 I I 

LU UJ UJ 

r* no o ** 

fn h o cr 

r~ o © cc 

cu cc p- © 

r- p* eg m 

<r m h ej 

fNJ p«- CO — « 

• • • • 

o o o o 


o ® O O 

■ I II 

UJ uj uj uj 

fw ® O' eg 

re r- qj rg 

p- a m eg 

cc ® o m 

h* f- eg p. 

m m -u 

rg r** **1 p- 

• • • • 

o o o o 


0 o o 

1 I I 

UJ UJ 111 

m p- r- o ® 
h r v g- qj 
g n ^ © cu 
eg m n* ,y m 

engj'Ugjxj- 

m oc r- m 
eg w* © «h 
• • • • • 

0 O o o o 

1 


o co <£ g <o 

v ^ co n- co 

O CC 1 IM «- <0 

g g p- r- fv 

co g) <4 n g 

xy in © © © 

o o o o o 


eg rg eg eg eg 

0 o o o o 

1 I I l I 
UJ LU UJ UJ LU 
a> in — < eg 
girgp- g 
»y in r- 

eg eg nj © in 
m m ^ O' 
rg tg rg h 

• # * • • 

o C o o c 


0 o 

1 I 


0 o 

1 I 


r- h y 
n- cr -< 
eg eg m 
r» r- >t 

CO h 

h- -o -• 


c 


• • 


c o 


m m & o e- — i cr 

•*T x#- rg o r*> cr -- 

r-xTr-rgo eg eg ro 

oc^iqo mn-xt 

h'CdjZ O x© O re 

•c n p. j • aip ^ 

m ci ^ — * x0 -« 

• • • cr rg • • • 

o O o o o C 


QOO 
I I I 
Uj uj uj 
® O eg 

in O' 

re co eg 

cc in m 
® eg cr 
m — < «© 
eg r^ in. 

• • • 

coo 


0 o o 

1 I l 

UJ UJ Uj 

p- »r eg ^ xr 

o m co r*» g? 

egaj—'cc^ 

o-nnee 

a; M *4 ^ IT 
<XJ^x0^4^ 
• • • • • 

0 o o c o 

1 


m © in p' p 

® it e- e* o 
eg o O' eg 
m pw eg in eg 

n- in •-< m ^ 
*r m ® <0 ® 

o o w o c 


rg rg eg rg eg 

0 o o o o 

1 I I I I 

LU LU Uj LU LL 

m m .-4 m a> 

o~ cj cd uj gj 

h in - y e 

C in nj a.' co 

tn xj- •> — if' 

pim-o? 

rg rg rj m »4 

• • • • • 

o o o c c 



o c o o 

II II 

UJ UJ LU UJ 

,$■ eg m eg cc m 

mar*- mop 

CD -* ® eg xT O' 

m >© e' rp px- r*. eg 

xj- in o ^ x»- m c 

eg — u im o j- 

^4 >o •— 1 in *-• r- uj 

• • • » • I c 

c o c o c o 


0 o o o c 

1 I ill 

LU UJ UJ UJ UJ 

x© rg m in xj- 

m o p- ao O' cc 

oh® px co r- 

xc g; m o — — 

xjin a in cr o 

eg —* in «-■ in 

-< o ~4 hioit 
• • • • • • 
O o o o o o 


- (\l -4 

0 o o 

1 I I 

U. UJ UJ 
IT. ® xj* rg 
in © r- o 
o e* >0 a* 
g eg «p M 

® O xj- xT 
rg rg m cj 

m 'O in — 
• • • • 

C o O C? 


rg 

xT 

xT 

xT 

m 

xj- 


o 


xj- ^ re eg o 
rg P- g m © 
g p h h gj 

y p ® < ^ 
in ,}• o ri g 

xj m *© xu >o 

-+ P-4 ** —• 

* • • • # 

C O O c » 


r\ r-4 

c o 
I I 
u- UJ 
cr xy rg 

o r* o 
n- gj cr 
rg xj- rg 

O xT xT 

rg «r. c 

xO u ' ■ 

CJ CJ c 


0 o 

1 I 

LU U 

CJ o o 

—• m xO 

in xj- cj 

O' o cj eg 
u- rv c; — 
O rg in H 
r h o il 
• • • a 
cj c cj 


eg _< 

?? 

UJ UJ 
oc xy eg 
•op V 1 
r- x© ® 
rg xy eg 
O xy xr 

eg m cj 
> c in 

O C CJ 


ej h >4 —• m 

o o o o o o o 

III I I I I 

LU Uj UJ LL UJ LU UJ 

O m m o®cc*^o 

-4 g h k o (*i h 

p in rn hop o p 

h o «o «o h rg v a 

m ^ *o m m — o 

rg x© o f^oo'tt.y 

xyxo<© m t— is ^ 

• •• 

O CJ O CJ O CJ CJ o 


cr a- ec a o 

eg g- «o h n 

«f o e eg h 

xy »x- eg ej 

xymomxy 

xy in ® »© x© 

* • ♦ • • 

o CJ CJ C CJ 


rg eg rg rg 

o o c o 

till 

UJ UJ Uj UJ 
H«®H 

pinm- 
a. in in rg 
xy m o u 

CL CL <© O' 
nrgHV' 
rg eg rg h 
• • • • 

CJ o c CJ 


eg eg rg eg nj 

C O O o c 

I l I » I 

u UJ UJ UJ UJ 

r- x© «— o ® 

P h gj y y 

o. © n h 
y CP xy -h 
eo xy o m m 

xy rg cj o 

MM (V IM H 
• • • • • 

CJ c o CJ c 


0 

1 

Uj 

•4- 

m 

oo 

*u 

•g- 

*g 


> 

cr m 


in 



2 


eg 


a 

o rg- 


eg 

0 

1 

UJ 

cc 

<© 

r- 

*v 

O' 

rg 

>© 

c 



eg fn 

0 o 

1 i 

LU UJ 

r- ac 

xy m 

e> p- 

n c 

O' xy 

r- m 

xy «o 

cj o 


m 

r- 

u. 

re 

CJ 


CJ 


LU 

>© 

IT' 

eg 

eg 

r- 

eg 

CJ 


eg 

O 


CJ 


pr, ^ rp — 

O C C C O 

III I | 

UJ UJ UJ LU U: 

xO O' in eg — 

a* o rg ^ m ® 

o r> c "xyu' 

— rg o- — ' C' X© 

m m ^h e> a xy rg 

u u cr rg r“ c. 

— • xy a — * -4 o a 

u r c c c' c/ 


> < 

<3 — — a 

rt C O a C 

a I l < l 

<3 LU U LL 

C C •—* «» — 

c -* r- -gu'p 
,-orgoo c eg r* 

cc-cuum— *rr,. 

o. c, eg o j xr ,%J n — 

2 up. fMPC-CC 
u xy o c — — 


coo 


e“. — — m ■ 

coo coo 

III III 
u: LU u> LU IU LL 

<D cr m r- in o 

cc. CJ eg — . ^ -4 

o p cj © y p 

— < rg O' o m *c 

p p h oc cr o- 

O cr iL eg xy 

— xy o «o x© o 


o c ® o 
i i l i 
UJ LL 1 Uj UJ 
in O C <© 

O X H P o> 
(MHlMOff 
h u o y a 
© © IV c xO 
eg m e-. «g m 
xy ,y xy o ~- 


c c u 


cce o o u c c 


d in o in O 
u n -4 cr eg 
ru p a o' in 
or a. xC o- cj 
fg tg cr rg -y 
^ in © •© g 


c* c o e c 


rg rv. fv rv rg 

0 o c; c o 

1 i < I i 

U 1 U- It LL u: 

a p ® p © 

© p p cr ® 

U' & P' u'. p- 

a «o n O eg 

e ® m g ® 

xy rv rg C O 

rv eg rv rg _4 

• « • • • 

o e c; c c 


a i i « i i t 

< U. u u U U, 

O c — xr xj- o 
C — rg* xy n cr 
rg c eg o n o it it- 

c. C- c 0-0 e - r - C 

y C rw O 1 'T C p- < 

s c p i^ p c p r : 

o x® u u ir> © 

c cu COO 


<llll cr 

a ui uj u. u.' < 

a © rg — •© eg cj 

C xfP "OP «r» 

•© eg — u <g c: a 

cj cr p a. P e - o. 

o xT cj a — • u- » 

5 e rg e ©. p g Q 

- to e - tr — — 

• ••••• 

CJ Cj g- c. CJ c 


< rg fv rv n rv e. 

a C C c- c O C 

a i i l I i l 

< L. LU LU U U- LL 

pevxymoc a c n o- g: a. 

<r >© xc n; r- o g © a p i 

r- orgxy a an <. © c p cv c: 

oc cr-or v g © L c. ©. c i 
— ‘ — C T ig r- X© U- rg c x* *© *© e 

xy in in o ® x© a if- eg rg o C- C 

rg rg ep rg — < 

OOCOCJCJ C, CJ C. c o c 


u Ll 


I < I 


I I 
u u 

c a. rv a © i\ a 

C C C r\ cr a t**. q <r c — 

o c Uj Li u u o. C u u U-O 

K UC-O - Cl-C P P I- -X 

g _'CPieut-CpaM\ 
C — 5 c Q r- Cl e* 


37 


7 BAC 75 E-C 2 C. 5 '.'I 6261 E -03 ' -0 . 3406031 F -03 - 0 . 27 3 6?2 2 F-C'i 


variables to be used in each subsequent iteration of the reduced-mass-flow solution. 
These maximum, minimum, and change values indicate to the user how 3(rV 0 )/3r and 

V 

F r are converging and give clues as to how to change DNEW and FNEW if convergence 
cannot be obtained. 

(5) This output corresponds to IMESH for one of the iterations of the reduced-mass- 
flow solution. The output is given along each horizontal mesh line at each mesh point. 
Typically, there are from 10 to 30 of these horizontal mesh lines, with 30 to 50 points 
in each. The output given at each mesh point includes the following: z- and r- 
coordinates of the point; stream function u; meridional velocity W m ; relative tangen- 
tial velocity W relative velocity W; critical velocity ratio W/W cr ; meridional flow 
angle a; relative flow angle /3; and mesh angle cp. 

(6) This output corresponds to ISLINE for one of the iterations of the reduced-mass- 
flow solution. This output is given along each streamline, corresponding to a given 
stream -function value. The points along the streamline correspond to where it is inter- 
sected by the vertical mesh lines. The origin for the m -coordinate is chosen so that 

m = 0 when z = 0. The output given at each streamline point includes the following: 
z-, r-, and m-coordinates of the point; meridional velocity W m ; relative tangential 
velocity W 0 ;. relative velocity W; critical velocity ratio W/W cr ; meridional flow 
angle a; relative flow angle /3; streamline curvature l/r c ; and within the blade -to- 
blade passage, the estimated suction-surface and pres sure -surface velocities. 

(7) This output corresponds to ISTATL for one of the iterations of the reduced- 
mass-flow solution. Instead of output being given along horizontal mesh lines (IMESH) 
or streamlines (ISLINE), it is now given in the other direction, along lines from hub to 
shroud. This output is given along each of the station lines specified by the input 
ZHST, ZTST arrays. Output is given at each point where these station lines are crossed 
by the streamlines (the values in the FLFR array). The output given at each point is 
identical to the output given for ISLINE, with the addition of the stream function u. 

(8) This output corresponds to ITSON. This is a listing of the information required 
to prepare input for the TSONIC blade-to-blade analysis program of reference 5. This 
information is printed in such a way that it corresponds to the input form for TSONIC, 
which is shown in figure 12. (The definitions of all the input variables for TSONIC are 
given in the TSONIC report, ref. 5. ) The TSONIC program was run for both the hub 
and tip blade sections of the compressor blade used in the numerical example of this 
report. The input numbers shown in figure 12 are those computed for use with the 
TSONIC run at the hub. 

Much of the information printed by MEREDL has to be corrected slightly before it 
can be used as input for TSONIC. However, the final three arrays - MR, RMSP, and 
BESP - can be used as they are printed, unless the user wishes to reduce the number 
of points in them. The user should recall that the MERIDL program takes into account 
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Figure 12. - TSONIC input form. Data shown are for numerical example at the hub. Inputs used for MR, RMS P, 
and BESP are the same as the data printed in item 8 of table 1. p. 33. 
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a loss in total pressure through the blade rows. TSONIC does not do this. Therefore, 
to ensure compatability between the programs, the BESP array calculated by MERIDL 
is reduced to reflect this loss in total pressure, both through the blade row and down- 
stream, before it is printed as input for TSONIC. 

The BESP array is calculated to correspond to the given mass flow, WTFL. Any 
value could be used for WTFL if BESP is chosen accordingly. The value calculated by 
MERIDL is 1 percent of the mass flow between two blades. 

Figures 13 and 14 indicate how the blade geometry information has to be corrected 
by the user before it can serve as input for TSONIC. First of all, the THSP1 and THSP2 
arrays will usually have to be corrected to the proper origin. As figure 13 indicates, 
the origin for the 0-coordinates is at the center of the leading-edge radius for TSONIC 
but at the mean camber line in MERIDL. 



Figure 13. - Blade leading edge showing origin of coordinates for MERIDL 
output and for TSONIC input. 


A layout (fig. 14) should be made to determine the remaining geometric input. To 
assist in making this layout, some additional coordinates are printed out. These addi- 
tional coordinates are tangential distances (r0) from a reference plane near the leading 
or trailing edge, against the m-coordinate. When these are plotted at the same scale, 
the leading- or trailing-edge radius can be plotted and measured. This will give RI1, 
RI2, ROl, and R02. Usually, RI1 is the same as RI2 and ROl is the same as R02. 

The layout, figure 14, indicates how the leading-edge radius should be measured. 

In addition, the tangent angles, BETI1 and BETI2, can be measured from the layout. 
Finally, the layout, figure 14, should be used to obtain the correction for the THSP1 
and THSP2 arrays. Figure 14 shows a tangential offset of 0. 00052 meter. This is con- 
verted to radians by dividing by the local radius of 0. 1537 to obtain a correction of 
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0. 00338 radians. Therefore, 0. 00338 was subtracted from all THSP1 and THSP2 co- 
ordinates. At the trailing edge, the 0-coordinate at the center of the trailing-edge 
radius is calculated from a similar layout. For this case, the trailing-edge 9 - 
coordinate is 0. 05886. Then 0: 00338 is subtracted to obtain 0. 05548. This is the cor- 
rect value to use for STGR. One further change should be made in the blade coordinates. 
The spacing between the MSP coordinates should be no closer than the original ZBL 
input for the blade shape. In this case, the ZBL coordinates were as much as 0.0011 
apart along the hub. Therefore, only every other point in the MSP and THSP arrays 
should be used. In other cases, only every third or fourth point may be used. Also, 
be sure that the first spline point after the tangent point is no clas er than the original 
ZBL input. Finally, be sure to change SPLNOl and SPLN02 when you do not use all the 
points. 

After all the geometrical input is obtained, all that remains is to choose the mesh- 
point numbers MBI, MBO, MM, and NBBI, as described in reference 5. The user 
may also wish to change some of the other input variables (ORF, REDFAC, DENTOL, 
and BLDAT to SURVL, for instance) since MERIDL only prints suggested values. 

The input used for the numerical example at the hub is shown in figure 12, except 
for the MR, RMSP, and BESP arrays. These three arrays were not changed from the 
output printed by MERIDL. 

When running the TSONIC program, minor changes to the blade geometry and tan- 
gent angles are often required to obtain a smooth blade surface and a desirable velocity 
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distribution. Changes of this type normally would not require the MERIDL program to 
be rerun. 

(9) This output consists of blocked and unblocked incidence and deviation angles at 
the blade leading and trailing edges where each of the horizontal mesh lines intersects 
the blade. The blocked incidence or deviation is based on the velocity diagram within 
the blade (subscript bf in fig. 15), whereas the unblocked is based on the free-stream 
(fs) velocity diagrams. This output is printed after each iteration in which output is 
given for either IMESH, ISLINE, or ISTATL. These angles are defined as shown in 
figure 15. For the blocked angles, the flow direction is corrected for blockage before 
the incidence and deviation are calculated. 



Figure 15. - Definition of incidence angles. 


(10) This output indicates the maximum relative change in relative velocity W at 
any point on the orthogonal mesh during an iteration of the reduced-mass-flow solution. 
When this value becomes less than the input value of VELTOL, the reduced-mass-flow 
solution is considered converged. 

(11) This output is analogous to outputs (5), (6), (7), and (8) for the variables IMESH, 
ISLINE, ISTATL, and ITSON but is given after the transonic velocity-gradient solution. 
Two solutions can be obtained by the velocity-gradient method, the larger or "super- 
sonic” and the smaller or "subsonic. " The input variable ISUPER controls which of 
these is obtained. The solution listed here is the smaller (i.e., the "subsonic" solu- 
tion). If REDFAC = 1.0, no velocity-gradient solution will be obtained or printed. 
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Plotted Output 


Since the printed output from a MERIDL run can be quite large, microfilm plots are 
made to enable the user to quickly check the quality of both his input and output data. 

The coding which generates these plots uses NASA Lewis in-house plot routines and 
would have to be recoded for operation on another facility. The principal sections of in- 
put data plotted are the upstream and downstream distributions of flow properties; the 
input blade sections; and the hub, shroud, and blade leading- and trailing-edge geometry. 
The amount of plotting of output data is controlled by the variable IPLOT. The plotted 
output data begin with the generated orthogonal mesh. Then for each outer iteration of 
the solution (indicated by IPLOT) and for the transonic solutions, streamlines are plotted, 
as well as mid-channel flow plane and blade surface velocities along each streamline 
from hub to shroud. The user should carefully check the plots of his input data to ensure 
that the program's spline fits of this data are smooth. 

Selected examples of some of the microfilm plots generated by the numerical ex- 
ample are presented in figure 16. The description of these plots follows: 

(1) Figures 16(a) to (e) present the input upstream and downstream flow conditions. 

These quantities (T!; p!; X or ; Pq or pressure loss; (rV 0 j or ) are 

all plotted against either stream function or radius, whichever was used as input. 

(2) Figures 16(f) and (g) indicate hub and shroud blade sections plotted from input 
data supplied by the user. All blade sections given as input are plotted. 

(3) Figure 16(h) shows the hub and shroud channel geometry and the blade leading 
and trailing edges. 

(4) Figure 16(i) shows the generated orthogonal mesh. 

(5) Figure 16(j) shows the streamline pattern for the final iteration of the reduced- 
mass -flow solution. 

(6) Figures 16(k) and ( l ) show meridional velocities and blade surface velocities for 
the hub and shroud blade sections for the same iteration. A similar plot is made along 
each of the streamlines from hub to shroud. 

Streamline and velocity plots similar to figures 16(j), (k), and (l) are repeated after 
each iteration of the finite-difference stream -function solution and also after the velocity- 
gradient solution. 
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(b) Inlet absolute total pressure. 

Figure 16. - Continued. 
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(d) Outlet absolute total pressure. 

Figure 16. - Continued. 
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(i) Orthogonal mesh in meridional plane. 
Figure 16. - Continued. 
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MERIDIONAL COORDINATE 

(k) Meridional and surface relative velocities. Streamline 1; normal 
ized stream function, 0. 

Figure 16. - Continued. 
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(i) Meridional and surface relative velocities. Streamline 11; normal 
ized stream function, 1.0000. 


Figure 16. - Concluded. 


Error Messages 


A number of error messages have been incorporated into the program. These error 
messages are listed here. Suggestions for finding and correcting the cause of the error 
are given. 


MM, MHT, NHUB, NTIP, NIN, NOUT, NBLPL, NPPP, NOSTAT, 

NSL, LSFR, LTPL, OR LAMVT IS TOO LARGE OR TOO SMALL 

The input dictionary gives the maximum and minimum value for all these input variables. 


LININT CANNOT FIND INTERPOLATED VALUE 
This message should occur only with erroneous geometry input. 


PASSAGE IS CHOKED AT THE INLET (OUTLET) WITH A MASS 
FLOW OF X. XXXX 

This message is printed if the choking mass flow at the upstream (downstream) bound- 
ary of the mesh region is less than the input mass flow (MSFL). Usually, the mass flow 
must be reduced, or there is an error in the input upstream (downstream) flow condi- 
tions. Also check hub and shroud coordinates at boundary. 


INLET (OUTLET) BOUNDARY CONDITIONS CANNOT BE OBTAINED 

This message is printed if the upstream (downstream) boundary conditions cannot be 
satisfied after 100 iterations. This may be caused by a severe variation of some up- 
stream (downstream) flow condition. Carefully check all input. 


PROGRAM STOPPED IN NEWRHO DUE TO EXCESSIVE STREAM FUNCTION 
GRADIENT 

This is usually caused by having supersonic flow in some region of the reduced- mass- 
flow solution. Try a smaller value for REDFAC. 
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UPSTREAM WHIRL OR TANGENTIAL VELOCITY IS TOO LARGE 


This message is printed if the upstream whirl gives a tangential velocity so large as to 
result in a negative relative stagnation temperature. The upstream whirl or tangential 
velocity given as input should be corrected. 


A VELOCITY GRADIENT SOLUTION CANNOT BE OBTAINED FOR VERTICAL 
ORTHOGONAL MESH LINE I = XX 

ANY SUBSEQUENT OUTPUT FOR THAT MESH LINE MAY BE IN ERROR 

If there is a problem in obtaining the approximate transonic solution for any particular 
mesh line, the attempt will be halted after 100 iterations, and this message will be 
printed. The problem may be caused by complex geometry or large gradients in up- 
stream or downstream flow conditions. Sometimes a finer mesh will help. 


MSFL EXCEEDS CHOKING MASS FLOW FOR VERTICAL ORTHOGONAL MESH 
LINE I = XX 

CHOKING MASS FLOW = X. XXXX 

This message is printed if the choking mass flow calculated by subroutine TVELCY is 
less than the input mass flow (MSFL). The following message will be printed after all 
vertical lines have been checked. 


CHOKING MASS FLOW IS LESS THAN THE INPUT MASS FLOW 

INPUT MASS FLOW = X. XXXX 

MINIMUM CHOKING MASS FLOW = X. XXXX 

A SOLUTION CAN BE OBTAINED ONLY IF INPUT MASS FLOW IS LESS 
THAN THIS MINIMUM CHOKING MASS FLOW 


If the previous message was printed for any vertical mesh line, then this message will 
be printed at the end to give the minimum choking mass flow, which is the choking mass 
flow for the blade row. 
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INRSCT HAS FAILED TO CONVERGE IN 20 ITERATIONS 
TOLERANCE = X. XXXX 

DISTANCE BETWEEN LAST TWO INTERSECTION POINTS = X. XXXX 

Subroutine INRSCT finds the intersection coordinates of mesh lines and streamlines with 
the blade leading or trailing edge by an iterative method. If the tolerance cannot be met 
after 20 iterations, the message is printed. If the distance between the last two inter- 
section points is only slightly larger than the tolerance, a satisfactory solution will be 
obtained, with some loss of accuracy. If the distance is excessive, there is probably 
some error in the geometry input. 


ROOT HAS FAILED TO LOCATE A ROOT IN THE INTERVAL (A, B) IN 
20 ITERATIONS 

ROOT ARGUMENTS -- A = X. XXXX B = X. XXXX 
Y = X. XXXX TOLERY = X. XXXX 

X FX DFX INF 

X.XXXX X.XXXX X.XXXX X 

i r i i 

Subroutine ROOT is called by subroutine MESHO to find the intersection between straight 
lines and spline curves in the process of generating the orthogonal mesh. If this mes- 
sage is printed, it is probably because the tolerance could not be met, although the 
problem may be caused by erroneous geometry input. 


NUMERICAL EXAMPLE 


An example of an axial-flow compressor rotor is used to illustrate the use of the 
program and to show the type of results which can be obtained. Flow is analyzed on the 
mid-channel flow surface of an axial-flow rotor designed with the computer program of 
reference 11. The design pressure ratio is 1.275, the inlet hub-tip radius ratio, 0. 5; 
the aspect ratio, 1. 5; the tip solidity, 1.0; and the tip relative Mach number at the inlet, 
0. 9. Although tip relative Mach number is near sonic, there are no locally supersonic 
regions on the meridional flow surface. 

The input for this example is shown in table II. The number of mesh points used 
was 861: 41 in the axial direction, and 21 in the radial direction. The printed output 
presented earlier is from this example. 
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TABLE n. - INPUT FOR AXIAL-FLOW COMPRESSOR ROTOR 
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TABLE Q. - Concluded. INPUT FOR AXIAL-FLOW COMPRESSOR ROTOR 
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REDFAC 

0.99999 

.70 




Finite difference (REOFAC =■ 0.999991 

Velocity gradient (REOFAC ■■ 0. 99999) 

Velocity gradient (REDFAC = 0.70) 


\ 


Blade leading edge — . 


Blade trailing edge- 


(a) Hub section. 

/Si 


V 

Blade leading edge-* 



Blade trailing edge 


tb) Mean section. 



Blade leading edge — 


Blade trailing edge 


z-coordinate, meters 
(c) Tip section. 

Figure 18. - Mid-channel velocities for axial-flow compressor example. 




This rotor was run twice, with values of reduction factor (REDFAC) equal to 0. 99999 
and 0. 70. This was done in order to test the effect of reduction factor and also the ac- 
curacy of the velocity-gradient method in comparison with the finite -difference method. 
Since both cases have REDFAC < 1. 0, both the subsonic (finite difference) and the 
transonic (velocity gradient) solutions will be obtained in each of them. (If 
REDFAC = 1. 0, only the finite-difference solution is obtained. ) 

The first case (REDFAC = 0. 99999) permits a comparison to be made between the 
approximate velocity-gradient method (used to obtain transonic solutions) and the more 
exact finite-difference stream-function solution (used for subsonic solutions) Since 
the reduction factor is 0. 99999 (<1. 0), both solutions will be obtained, but the answers 
should be very close if the methods compare well. 

The second case (REDFAC = 0.70) permits a comparison of reduction-factor effect. 
Since the finite-difference solution is now obtained at only 70 percent of the mass flow 
and rotational speed (see appendix E), the results for the full-mass-flow velocity- 
gradient solution will be less accurate than when a higher reduction factor was used. 

Figure 17 shows the streamline plots for the final finite-difference iteration for both 
the REDFAC = 0. 99999 and REDFAC = 0. 70 solutions. The differences in streamline 
location are small, but program printout of streamline curvatures shows that these cur- 
vatures differ significantly. This illustrates the point that as high a value of REDFAC 
should be used as possible. 

This fact is further illustrated in figure 18, which shows mid-channel flow surface 
velocities for both reduction factors and both types of solution, the finite difference and 
the velocity gradient. On all parts of this figure, the solid line represents velocities 
from the final finite-difference iteration of the REDFAC = 0. 99999 solution. This solu- 
tion is mathematically most accurate. 

Figure 19 shows a comparison of blade surface velocities estimated from the 
MERIDL program (see appendix G) and calculated from the TSONIC program (ref. 5). 

The TSONIC velocities are more accurate for the design blade shape. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 16, 1972, 

501-24. 
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APPENDIX A 


GOVERNING EQUATIONS 


The coordinate system is shown in figure 2. Since 0 is a function of r and z on 
a meridional stream surface, 9 can be eliminated so that the two independent variables 
are r and z. 

For the subsonic solution, the stream function is used. The stream function u 
used herein is related to the stream function defined in reference 1 by i// = uw. 

With this substitution in equation (107a) of reference 1, we obtain in appendix B the 
basic differential equation which must be satisfied by the stream function under the 
given assumptions: 


3 2 u ^ 3 2 u 
0Z 2 3r 2 


i 

r 3r 


I 1 VB + 1 Vp\ • Vu + rBp 

W 0 3(rV 0 ) 2 1 

y + r + f^ 

\B p / wW z 

1 

QJ 

H 

I 


= 0 


(Al) 


where 




1 , 

(r 

3p” 

_ J_ 

31 

2C p 

^P" 

3r 

fprr 

3r 



(A2) 


. 2 
$ = w r 


RT” 8p M 
p” 3r 


(A3) 


F = — - ^2 (A4) 

r dr p 39 

Equations (Al) to (A4) are derived from Wu's equation in appendix B. Note that all 
the partials in equations (Al) to (A4) are on the stream surface except for dp/39 in 
equation (A4), which is at constant z and r. 

The derivatives of the stream function satisfy the equations 


3u rB P W r 

3z w 


(A5) 


3u _ rB P W z 
3r w 


(A6) 
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For the final transonic solution the velocity-gradient method is used. The stream- 
line curvatures and flow angles needed for the velocity-gradient equation are obtained 
first from a reduced -flow subsonic -stream -function solution. 

The velocity-gradient equation is 


dW = 


aW + b + — + d cos /3 )dt + — + Wf 


W 


W 


(A7) 


The coefficients, a, b, c, and d, are given by different experessions in the blade region, 
in the upstream region, and in the downstream region. These coefficients are given as 
follows: 


Blade-region coefficients: 


= c°g!g c oslg; _ sin 2 g c_ 0 _s_^ + sin a sin p cos pM 
r r 3t 


dW /dW 0 i 30 i 

b = cos 0 — sinitv - w) - 2u> sin 0 cos <p + r cos 01 + 2w sin a) — > (A8) 

dm ydm /at 

c = 0 


d = 0 


Upstream-region coefficients: 


a 


_ cos (a - <p) 
r 


b =0 


(x - u>r 2 \ 

2 2 
— - ^ cos(cy - (p) + — + a>r - cos <p 

\ r 2 ) 

r c r 


dW 

d = si^cy - cp) 

dm 


y (A9) 
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Downstream-region coefficients: 


_ coslo 1 - cp) 
r 


c = - 



K)„ - 


b = 0 
ur 2 


cos(cy - (p) + 


K), 


+ cor 




(A10) 


cos cp\ 


dw 

d = sin(cy - cp) 

dm 


Finally, in all three regions, we have 

DOM t 

e = C p T| - co dX + CpdT" - £±-dp" 

f - Rd P" _ dT M 
2C p p" 2T" 

Equations (A7) to (All) are derived in appendix C. Equation (A7) is solved as an 
initial-value problem, where the initial value of W is specified at the hub for any given 
t-line running from hub to tip. By finding several solutions for varying values of W at 
the hub, a solution satisfying continuity will be found; that is, the solution will satisfy . 

r t 

pW rB cos ( a - cp) cos 0 dr = w (A12) 

When equation (A7) has been solved, subject to satisfying equation (A12), for every hub- 
to-tip mesh line in the region, the entire velocity distribution is obtained. 

The solution obtained by either the finite-difference or velocity-gradient method is 
for the mid-channel surface between the blades. Of greater interest are the blade sur- 
face velocities. These can be estimated since the blade loading depends on the rate of 
change of whirl. By assuming a linear variation of velocity between blade surfaces, 



-\ 


(All) 
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velocities can be calculated, 
surface velocities is derived: 


In appendix G the following equation for calculating blade 


R d(rVJ 

W, = W . . - ^ cos jS - 

1 mid 2 dm 


R ' d < r V 

W tr = W mld+ |cos^-_l 


J 


(A13) 
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APPENDIX B 


DERIVATION OF STREAM-FUNCTION EQUATION 


Wu derives the following equation (eq. (107a), ref. 1) for the stream function on a 
meridional stream surface: 


3z 2 3r 2 r 3r \B p / dty/dr 


W 0 3 < r V 31 , T 3s , F 


= 0 


3r 


3r 


3r 


(n 


The partial derivatives here all refer to the rate of change on the meridional stream 
surface. Wu used a bold partial derivative sign to indicate this, but it is not necessary 
here since we consider the equation only on the stream surface. The quantity B is 
proportional to the blade -to-blade stream surface thickness. Since the variation in 
local stream sheet thickness is not known, the overall blade-to-blade spacing is used. 
The quantity B is a function or z and r, so that we can use 


B(z, r) = 0 tr (z, r) - 9 t (z, r) 


(B2) 


where 0 tr and 9^ are the 9 -coordinates of the trailing and leading blade surfaces. 

Equation (Bl) is a nonlinear partial differential equation. The solution can be ob- 
tained by solving a linearized form and then making corrections to the nonlinear terms 
to improve the solution. After several iterations, the true solution to the nonlinear 
equation is obtained. The equation can be put in a linear form by expressing dip /dr 
in the denominator in terms of W and using velocities and densities from the previous 
iteration. For the first iteration,, the nonlinear terms are omitted. 

First, the stream function will be normalized to be 0 at the hub and 1 at the tip. 
This is done by letting 


i|/ = uw 

where u is the normalized. stream function. Then 

3u _ rB P w r 
3z w 

3u . rB ? W z 
3r w 


(B3) 


(B4) 


(B5) 
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By using equations (B3) and (B5), equation (Bl) can be written 


a 2 u t d 2 u 


3z 3r‘ 


1 3u 
r 3r 


(L VB + - Vp\ 
\B P / 


Vu + 

wW z 


r 3r 


- — + T — + F 
3r 3r 


= 0 

(B6) 


Rothalpy I is defined by 


I = CpT| - cuX 


(B7) 


The entropy change ds can be calculated from 


ds = C„41 .RdE 

p T p 


Since the entropy of a particle at actual flow velocity is the same as the entropy of 
the particle at stagnation conditions (either absolute or relative), we can also use 


ds = C. 


dT” 


R dp" 


(B8) 


This takes less calculation than using static values. 

Equation (97) of reference 1 gives F r (the radial component of the vector F) as 


n^r p 30 


(B9) 


Note that 3p/30 is at constant r and z, and not on the stream surface as the other 
partials are. The unit normal vector n to the stream surface is normal to any curve 
on the stream surface, so that 


n r dr + n^ r d0 + n z dz = 0 (BIO) 

In particular, at constant z, we have 


n 


r 


n 


0 



(Bll) 
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where r(dd/dr) is the tangent of the blade lean angle at constant z. Using equation (Bll) 
in equation (B9), we get 


F 38 1 dp 
r dr p 39 


(B12) 


The static temperature T is calculated from the velocity W by 


T _ T , W 2 + 2o)X - (o;r) 2 
‘ 1 2C p 


(B13) 


When W is zero, equation (B13) gives the relative stagnation temperature, T". 

Equations (Al) to (A4) are obtained from equations (B6), (B7), (B8), (B12), and 
(B13). 


STREAM-FUNCTION EQUATION IN s- AND t-COORDINATES 

The solution of equation (Al) is obtained by using finite-difference methods on an 
orthogonal mesh. The orthogonal mesh coordinates are s in the throughflow direction 
and t in the hub-to-tip direction. The derivatives of the stream function in the s and 
t directions are related to the velocity components by 


w 

3s w 


= w 

3t w S 



(B14) 


Since s and t coordinates are used in solving equation (Al), the partials must be 
expressed as partials with respect to s and t. This is done by using the chain rule: 

_L _ jLi§ + ±dt 

dr ds dr 3t 3r 


3 3 

= sin <p — + cos <p — 

3s at 

and 

3 3.3 

— = cos cp sin cp — 

3z 3s 3t 
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Now, apply this to the second partial of u to obtain 


a 2 u t a 2 u _ a 2 u ; 3 2 u t au a cp au d<p 

dz 2 dr 2 as 2 at 2 as at at 0S 


(B15) 


In the program sin cp and cos cp are stored, so that it is more convenient to calculate 
the partials of cp by 


dcp _ 

l 

3(sin cp) 

as 

cos cp 

ds 

dcp _ 

1 

3(sin cp) 

at 

cos cp 

at 


(B16) 


The s, t orthogonal coordinate system is in the same length units as z and r. In this 
case, the gradient can be expressed directly in s and t coordinates, that is, 


VB • 


3B 9u 9B 0u 
ds ds dt dt 


(B17) 


Finally, putting all this together, we have 
" 2 u + 9 2 u _ 3u fsin cp + 1_ 3B + 1 d£ _ 1 8(sin (p)" l 

3s 2 at 2 9s L r b as p ds cos cp at J 


aufcos cp 

1_ 3B 

, 1 S£ + 1 

3(sin cp) 

atL 

r 

b at 

p at cos cp 

ds 

rBp j 


sin cp 

3(rV 0 ) 

I- pon cn 

9(rV 0 ) 

wW z | 

r 

ds 

at 


+ tw* + ? + F r = 0 


(B18) 


where and F r are as defined in equations (A2), (A3), and (A4) 


CALCULATING F r AND PARTIALS OF (rV 0 ) 

In equation (B18), F r and the partials of r Vq must be calculated differently out- 
side the blade and within the blade. Within the blade, the partials of rVg are calculated 
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by a finite -difference approximation using the previous iteration. For the first iteration, 
the partials are assumed to be zero. Outside the blade, the partials are calculated from 
the input whirl distribution. From the input, the whirl is specified (or can be estimated) 
as a function of the stream function u. Hence, upstream of the blade, 

8 ( r V _ 3\ _ d\ 3u = dA rB ^ W z ( B19 j 

3r 3r du 3r du w 


Similarly, downstream of the blade, 

9(rV g ) _ d ( rV g) Q rBpW z 
3r du w 


(B20) 


The F f is caused by the radial pressure gradient induced by the blade lean. We 
can calculate F f from equation (B 12). The blade-to-blade pressure gradient 3p/30 
can be calculated from the blade loading by assuming constant entropy blade to blade, 
so that 


dp = pC dT (B21) 

r 

From equation (B13) we get 

iE = - pW — (B22) 

3 9 30 


since T!, A, and r are all constant from blade to blade. 
Now, substitute this in equation (B12) to get 

' F_ = - — W M 

. 3r 30 


(B23) 


The blade-to-blade velocity gradient is calculated from dW^/dm, as explained in ap- 
pendix G, by using equation (G2). Outside the blade, there is no blade-to-blade pres- 
sure gradient. Hence, 


F 


r 


= 0 


(B24) 


outside the blade. 
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APPENDIX C 


DERIVATION OF VELOCITY-GRADIENT EQUATIONS 


The general velocity-gradient equation is an expression for the value of the direc- 
tional derivative of the relative velocity. The velocity-gradient equation is derived 
from Newton's three-dimensional force equation. In this program the velocity-gradient 
equation is used in two ways. One way is to obtain the flow distribution at the upstream 
and downstream boundaries to establish the boundary conditions for equation (Al). The 
second way is to determine the approximate velocity distribution when there is locally 
supersonic velocity on the meridional mid-channel stream surface. 


VELOCITY -GRADIENT EQUATION FOR DETERMINING BOUNDARY CONDITIONS 


The general velocity-gradient equation can be written 


where 


™ = a 4l + b^ + c^ + A(C -i- M ®: + T^ 

dq dq dq dq W \ p dq dq dq. 


(Cl) 


a - ^ co -£- a c . ? s2 g - + sin a cos 0^51 - 2n> sin 0 

r r dm 


>1 


? dW 

u W cos p sin a . „„ _ _ „ m 
b = - + cos at cos p 

r„ dm 


(C2) 


m Q 

c = W sin a sin p cos p + r cos P[ + 2u> sin «) 

dm 


This is the same as equation (B13) and (B14) of reference 2, but with the additional 
term, T(ds/dq), to allow for variation in entropy. 

Outside the blade row, c = 0. This follows from the fact that d(rVg)/dm = 0 in 
this region. Using this identity and the relation Wg = Vg - cur, we can derive the ex- 
pression 


dWg Wg + 2 cur 

sin a 

dm r 
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When this is substituted in the expression for c, we find that c = 0. 

It is reasonable to assume that the boundary is normal to the flow along the upstream 
and downstream boundaries. Then q = n and 


— = cos a 
dn 


> 


U4i 

— = - sin a 

dn J 


(C3) 


When this is substituted in equations (Cl) and (C2), we obtain 

g . W cos 2 g . W slnWj* . 2m cos „ stn p.l / c ^. M to +T ds\ (c4) 
dn r c r W l p dn dn dn/ 

This equation is the basic velocity-gradient equation along a meridional streamline nor- 
mal. At the upstream boundary (or any place upstream of the blade) X, T!, and p| are 
specified as a function of the stream function. Also, along the boundary, a is assumed 
to be known, since it is assumed that the boundary is normal to the streamlines. The 
angle /3, however, is not known directly. Therefore, in equation (C4), we desire to 
eliminate in favor of X, r, and W. This can easily be done since 



- cur 


(C5) 


and 


W B 

sin f3 = — 
W 


(C6) 


Then cos 0 is obtained from 

cos 2 p = 1 - sin 2 /? (C7) 

The entropy change, ds/dn, can be calculated from the differential form of the 
second law of thermodynamics: 


ds _ C p ~ R dT _ R dp 
dn T dn p dn 
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Since the entropy of a particle at actual flow velocity is the same as the entropy at stag- 
nation conditions, we can also use 


ds_(V R ) dT i R d Pj 

dn T! dn p[ dn 


(C8; 


The meridional-plane streamline curvature l/r c is assumed to vary linearly from hub 
to tip. 

By substituting equations (C5) to (C8) in equation (C4) we obtain the velocity-gradient 
equation in the form used in the computer program for the upstream boundary: 



(C9) 


where 


a = 


b = - 


[A - cur 


X - cur cos a /x 2x 
(A + cur ) 


c = C dT! - cu dA 
a p i 


RT 1 

% - < c p - R > dT ! - -T dp i 


(CIO) 


Similarly, we can obtain a velocity-gradient equation for the downstream boundary. 
Here we use downstream conditions for temperature, density, and whirl, so that the 
downstream velocity-gradient equation is equation (C9) with coefficients given by 
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RT 1 

Cb = (Cp - R)dT’ - -2 dp’ 

Pn 


VELOCITY-GRADIENT EQUATION FOR CALCULATING 
APPROXIAAATE TRANSONIC VELOCITIES 


We start with the general velocity-gradient equations (Cl) and (C2). These equa- 
tions will be applied along vertical mesh lines in the t-direction, so that q = t. (See 
appendix B for a description of the orthogonal mesh. ) Then 


— = cos cp 
dq 




> 


Uii . i 

— = - sin cp 

dq J 


(C12) 


Note that, by using equation (C12), 


■ / \ - dr dz 

sin(of - cp) - sm a — + cos a — 


dq 


dq 


cos(a - cp) = cos a — - sin a — 


dq 


dq 


(Cl 3) 


Using equations (C12) and (Cl 3) in equations (Cl) and (C2), we can get 
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dW „ Wcos 2 ffcos(c»-y) + cos p sln(a _ v) . (wsir?0 + 2u) sln ^L, ? + c d£ 
dt r c dm \ r dt 


+ — (C n — - co — + T— I 
W \ p dt dt dt 


(C14) 


We can express T ds in terms of the relative velocity W and relative stagnation 
conditions. We use 


By using equations (B8) and (C15), 
T ds = C p dT" 


T = T” 


we get 


RT" 

P" 


dp” 


w 2 


2C P 


W 2 / 

'r dp" 


^2p" 



Now we can write the velocity-gradient equation (C14) as 


dW = (aW + b) dt + — + Wf 
W 


(Cl 5) 


(C16) 


(C17) 


where 


= cos 2 fl cos (a - cp) _ sin 2 g_cqs £ + sin a sin p CO s 0 ^ 


dt 


dW„ 


/dW „ 


’ m l ’ 7 0 \ d0 

b = cos /3 — sin ( a - cp) - 2co sin /3 cos cp + r cos 0 1 + 2co sin c* — 

dm \ dm / dt 


» _ RT" „„„ 


y (CIS) 


e = C dT! - co dA + C dT” - dp 
PI P „»» 


f _ R dp" _ dT^ 
2C p p” 2T” 




Equations (C17) and (Cl 8) can be used directly within the blade row. However, 
outside the blade row, the angle j3 is not known directly, so that we desire to eliminate 
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/ 3 in favor of X, r, and W, as was done on the upstream boundary. For this, we use 
equations (C5) to (C7). Also, c = 0 in equation (C2) outside the blade as was shown 
previously. With these substitutions, equation (C14) becomes 


dW = 


a w + — + d cos p dt + — + Wf 


W 


W 


(C19) 


where 


a _ cos (a - (p) 


c = - 


fX - cor 


2 2 

X - cor / \ X + cor 

cos( a - (f) + cos (p 


dW 

d = sin(cy - cp) 

dm 


e = C dT| - co dX + C dT ,f - dp” 
Pi P n" 


f = - dT” [ R dp” 
2T" 2CpP” 


(C20) 


cos /3 = 'l 1 - 


X - cor 
i rW 


j 


Downstream of the blade, outlet whirl and absolute stagnation temperature and pres- 
sure should be used. Hence, in equation (C20), is used instead of T!, p' o is used 
instead of p!, and ^rV^j is used instead of X. Equations (A7) to (All) are obtained 

directly from equations (C17) to (C20). 
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APPENDIX D 


LOSS CORRECTIONS 


An approximate loss correction is made by reducing the relative stagnation pressure 
used in the program. This loss in stagnation pressure is specified by the input either 
directly or by specifying outlet stagnation pressure> These quantities may vary from 
hub to tip. If the outlet stagnation pressure is given as input, the fractional loss of 
stagnation pressure is calculated. So, in all cases, the fractional loss of stagnation 
pressure at the outlet is known from hub to tip. This loss is assumed to be zero up- 
stream of the blade, then to vary linearly from leading edge to trailing edge within the 
blade, and to be constant downstream of the blade. 

For the case where the downstream stagnation pressure is given as input, the frac- 
tional loss of stagnation pressure is calculated. First, the outlet total temperature T^ 
is calculated from the change in whirl along a streamline, 


T o 




(Dl) 


Then the ratio of actual to ideal stagnation pressure is calculated from 


p; 



ideal 



(D2) 


Within the blade, the loss is distributed linearly from zero at the leading edge to the 
fraction given or calculated at the blade trailing edge. The loss fraction is the same 
whether it is expressed in terms of relative or absolute stagnation pressures. That is, 

-jell- = - 2 l_ (D3) 

p ideal p ideal 


When this ratio is known, the density can be calculated, since 


P 




p ideal 


(D4) 
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Equation (D4) is derived as equation (B6) in reference 13. The value of P"/p-^ ea q varies 
at each point of the region but does not change after the initial calculation. The value 
of p is used in checking continuity by means of equation (A12). Therefore, in the pro- 
gram, the loss correction is made by reducing the tangential blade space B at each 


mesh point. 
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APPENDIX E 


DEFINING REDUCED-MASS-FLOW PROBLEM 


When the mid-channel meridional-plane solution has locally supersonic flow, the 
solution cannot be obtained directly by solving the stream- function equation (eq. (Al)). 
However, an approximate solution can be obtained by getting a reduced-flow solution and 
extending this to the full-flow solution by the velocity- gradient method. This full-flow 
solution depends strongly on the reduced-flow solution, so that it is important to estab- 
lish the conditions which will give the most suitable stream- function solution. With a 
reduced mass flow, if both the inlet and outlet whirl and the rotational speed are reduced 
in the same ratio as the mass flow, a similar flow will result. That is, the reduced- 
flow streamlines will have approximately the same angles and curvatures as the full- 
flow streamlines throughout the flow field. 

Another consideration is the boundary conditions at the upstream and downstream 
boundaries. Again, the flow distribution will be similar to the reduced mass flow if the 
mass flow, inlet and outlet whirl, and rotational speed are all reduced by the same 
ratio. 

In summary, then, for the reduced-mass-flow solution, the mass flow, inlet and 
outlet whirl, and rotational speed are all reduced in the same ratio. This reduction 
factor is specified in the input by RED FAC, as explained in the section INPUT. 
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APPENDIX F 


INCIDENCE AND DEVIATION CORRECTIONS 

The solution region is divided into three subregions: upstream, within the blade, 
and downstream. These three regions should match each other, that is, there should 
be no difference in the tangential momentum of a particle on either side of the lines at 
the leading and trailing edges. If there were no incidence or deviation, this condition 
would be satisfied automatically by using the actual blade shape in the blade region. 

Since there is almost always incidence and deviation, the mean flow will not follow the 
blade surface near the leading and trailing edges. Because of this, an empirical correc- 
tion is made to the blade shape near the leading and trailing edges. 

The blade shape correction is made for a distance from the leading or trailing edge. 
This distance varies between 1/6 and 1/2 of the blade chord, depending on solidity. 

The solidity is defined to be the ratio of the true blade chord along a streamline to the 
blade-to-blade distance. When there is a change in radius, the solidity differs at lead- 
ing and trailing edges. For calculation purposes in the program the true blade chord 
along an s-coordinate line is used instead of that along a streamline. For a solidity of 2 
or more, the distance is 1/6 of the chord; and for a solidity of 1/2 or less, the distance 
is 1/2 of the blade chord. The distance varies linearly with solidity between these 
limits. The correction is made to the blade angle so that tangential momentum is con- 
tinuous at the blade leading or trailing edge, and then the angle correction is linearly 
decreased over the prescribed distance within the blade. 

The free- stream flow angle is not known in advance at the blade leading edge. 
Therefore, an iterative procedure is used to make the blade shape correction. On the 
first iteration, no blade shape correction is made. After the first iteration, an approxi- 
mate solution is obtained, and the flow angle /3 is calculated throughout the region. 

From the requirement of continuous tangential momentum, 

w fs ■ w M 

From continuity, 

(w ) -h- = (wb) 

' m 'fs NBL ' m 'bf 


Hence, since tan j3 = W 0 /W m 
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(FI) 


tan = tan /9 fg 


/ B x NBL \ 
\ 2tt / 


It is assumed that the s-coordinate is close to a streamline. Therefore, only dd/ds 
will be changed, and 39/ dt will not be changed. We can express (3 in terms of the 
partials of 9 by 

tan /3. t = r = r ^ cos^ - (p) + — sin(® - cp)\ 
dm \\ds/ hf 9t -I 

When this is solved for (30/3s) b f, we obtain 



tan /3jj f 
r 


- — sin(£y - <p) 

at 


cos(a - <p) 


(F2) 


Equations (FI) and (F2) are used to calculate (30/3s)^ at the blade leading and 
trailing edges. The difference (30/3s) bf - (30/3s) b is then varied linearly within the 
blade for the specified distance. The calculations for these blade shape corrections are 
done by subroutine INDEV. 
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APPENDIX G 


BLADE SURFACE VELOCITIES AND BLADE-TO-BLADE AVERAGE DENSITIES 


The blade surface velocities can be calculated once a meridional mid-surface solu- 
tion is obtained. The blade -to -blade velocity gradient 9W/30 depends on dW 0 /dm. To 
obtain the desired relation, we use the velocity-gradient equations (Cl) and (C2). In 
this case, q = 9 and dr/d0 = dz/dd = 0. Also, it is assumed that Tj, X, and the en- 
tropy s are constant in the blade-to-blade direction. With this, equations (Cl) and (C2) 
become 


3W 

= W sin a sin /3 cos /3 + r cos /3 

89 



+ 2oj sin cv| 


(Gl) 


Using the fact that V 0 = W 0 + cor and W g = W sin 0, we can rewrite equation (Gl) as 


aw 


= cos 
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d(rV fl ) 

0 — - 

dm 


(G2) 


Since only a mid-channel solution is obtained, we assume that cos 0 and d(rV 0 )/dm 
are constant from blade to blade. This means that W varies linearly from blade to 
blade. (This may be in considerable error near the leading or trailing edge. ) With this 
assumption, then 


W mid = 


+ W tr 


(G3) 


Integrating (G2) from blade to blade (from 9 = 9^ to 9 = 0 t ) and using (G3), we obtain 


W tr = W mid 


W * = W mid 


. B cos j3 d ( rV e) 
. , + 


2 dm 

B cos d(r V 
2 dm 


> 


(G4) 


Note that equation (G4) is very close to that developed by Stanitz (eqs. (16) and (17), 
ref. 16). 
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Equation (G4) is used in subroutine BLDVEL to obtain blade surface velocities from 
the mid-channel meridional solution. 

It is desirable to consider the blade-to-blade variation in density to satisfy continu- 
ity for the blade passage. Equation (G4) gives information we need to do this. Since the 
mid-channel solution is considered to be representative of all meridional-plane stream 
surfaces, we can consider the solution to be based on average blade-to-blade conditions. 
This means that in equation (B14) average blade-to-blade velocity and density should be 
used. Equation (B14) becomes 


— - - — (pw t ) 

3s w v 'av 




-=-K) 

at w ' SA 


(G5) 


av J 


If we use Simpson's rule to calculate the average value, we obtain 


/ W ) = p l W s, l * 4p mid W s. mid + p tr W s, tr 

' s 4v 6 


(G6) 


The velocity component W„ can be expressed in terms of the velocity and flow angles by 

S 


W = W cos /3 cos(o' - i p) 

b 


(G7) 


This relation holds also on the blade surfaces, since the flow angles are assumed to be 
constant in the blade-to-blade (6) direction, so that 


W s i = Wj cos /3 cos(a - tp ) 

W gj tr = W tr cos 13 cos (at - (p) 

By using equations (G3), (G4), and (G8) in equation (G6), we can express as 


(G8) 


( pW s) = p av W s + ~ — — cos P cos( " • ' W tr } 

3.V 12 


(G9) 


where the subscript mid is omitted from W g , and 
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p 


av 


p l + ^mid + p tr 
6 


When a solution of equation (Al) is obtained, all quantities in equation (G9) can be esti- 
mated from the previous iteration, except for W . So we solve for W in equation (G9) 

s s 

and use equation (G5) to obtain 


W 


s 


at ( P L - P tr )cos ]3 cos(a - (p)(W t - W tr ) 

rB ^av 12 ^av 


(G10) 


This is the equation used in subroutine NEWRHO to calculate W g from 3u/3t. On the 

first iteration, p_ r = p! is used, and the second term is omitted. After this, the val- 
a,v 1 

ues of p mi( j, P tr > Pi, W i, W tr , cos /3, and cos(cy - <p) from the previous iteration are 
used. 

In a parallel manner, we can derive the equation for W^.: 


w 


3u 


w t = - 


3s ( p l " PtP C0S & sin ^ Q ' " VK W l " W tP 


rBp 


av 


12p 


av 


(Gil) 
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APPENDIX H 


SYMBOLS 

a coefficient in velocity-gradient equations (A7) and (Cl) 

B tangential space between blades, rad 
b coefficient in velocity-gradient equations (A7) and (Cl) 

C specific heat at constant pressure, J/(kg)(K) 

r 

c coefficient in velocity-gradient equations (A7) and (Cl) 
d coefficient in velocity-gradient equation (A7) 

e coefficient in velocity-gradient equation (A7) 

F vector normal to mid-channel stream surface and proportional to tangential pres- 
sure gradient, N/kg 

f coefficient in velocity-gradient equation (A7) 

2 2 

I rothalpy, CpTj - cuA., meters /sec 

m meridional streamline distance, meters 

n unit vector normal to mid-channel stream surface 

2 

p pressure, N/meters 

q distance along an arbitrary space curve, meters 

R gas constant, j/(kg)(K) 

r radius from axis of rotation, meters 

r radius of curvature of meridional streamline, meters 
c 

S entropy, J/(kg)(K) 

s distance along orthogonal mesh lines in throughflow direction, meters 
T temperature, K 

t distance along orthogonal mesh lines in direction across flow 
u normalized stream function 

V absolute fluid velocity, meters/sec 

W fluid velocity relative to blade, meters/ sec 

w mass flow, kg/sec 

z axial coordinate, meters 
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a angle between meridional streamline and axis of rotation, rad; see fig. 2 
/3 angle between relative velocity vector and meridional plane, rad; see fig. 2 
y specific -heat ratio 

£ coefficient in stream -function equation, defined in eq. (A3) 

9 relative angular coordinate, rad; see fig. 2 
A prerotation, ^rV^ , meter s^/sec 

coefficient in stream -function equation, defined in eq. (A2) 

3 

p density, kg/meter 

cp angle between s-coordinate line and axis of rotation, rad; see fig. 4 
\j/ stream function, kg/sec 

w rotational speed, rad/sec; see fig. 2 

Subscripts: 

av average blade-to-blade value 

b blade 

bf blade flow 

cr critical 

fs free stream 

h hub 

i inlet 

l blade surface facing direction of positive rotation 
m component in direction of meridional streamline 

mid mid-channel blade to blade 

o outlet 

r component in radial direction 

s component in the s-direction 

t tip, or component in the t-direction 

tr blade surface facing direction of negative rotation 
z component in axial direction 

9 component in tangential direction 
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Superscripts: 

' absolute stagnation condition 

" relative stagnation condition 
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